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ABSTRACT 

We present the results of a near-IR spectroscopic survey of 110 Class I proto- 
stars observed from 0.80 fim to 2.43 /mi at a spectroscopic resolution of R=1200. 
This survey is unique in its selection of targets from the whole sky, its sample 
size, wavelength coverage, depth, and sample selection. We find that Class I 
objects exhibit a wide range of lines and the continuum spectroscopic features. 
85% of Class I protostars exhibit features indicative of mass accretion, and we 
found that the veiling excess, CO emission, and Br 7 emission are closely related. 
We modeled the spectra to estimate the veiling excess (r fc ) and extinction to each 
target. We also used near-IR colors and emission line ratios, when available, to 
also estimate extinction. In the course of this survey, we observed the spectra 
of 10 FU Orionis-like objects, including 2 new ones, as well as 3 Herbig Ae type 
stars among our Class I YSOs. We used photospheric absorption lines, when 
available, to estimate the spectral type of each target. Although most targets 
are late type stars, there are several A and F-type stars in our sample. Notably, 
we found no A or F class stars in the Taurus-Auriga or Perseus star forming 
regions. There are several cases where the observed CO and/or water absorption 
bands are deeper than expected from the photospheric spectral type. We find 
a correlation between the appearance of the reflection nebula, which traces the 
distribution of material on very large scales, and the near-IR spectrum, which 
probes smaller scales. All of the FU Orionis-like objects are associated with re- 
flection nebulae. The spectra of the components of spatially resolved protostellar 
binaries tend to be very similar. In particular both components tend to have sim- 
ilar veiling and H2 emission, inconsistent with random selection from the sample 
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as a whole. There is a strong correlation between [Fe II] and H2 emission, sup- 
porting previous results showing that H 2 emission in the spectra of young stars 
is usually shock excited by stellar winds. 

Subject headings: surveys, stars:formation, stars: pre-main sequence, infrared: 
stars, techniques: spectroscopic 



Introduction 



Pre-main sequence stars have been classified by the slope of their spectral energy dis- 
tributions (Adams et al. 1987, Lada 1990), which traces the temperature and distribution 
of circumstellar material around the young star. Among these, the T Tauri (Class II) stars 
are the best characterized. They are optically visible, allowing many to have been discov- 
ered and spectroscopically studied since they were first identified as pre-main sequence stars 
(e.g. Walker 1956, Herbig 1962). The Class I objects are characterized by having increasing 
flux density with wavelength in the infrared due to their circumstellar envelope and disk, 
and are optically obscured. Comparatively few Class I objects are known, and they have 
only been st u died in detail with the development of infrared array detectors. Most recently, 
Evans et al.l ( 120091 ) characterized the YSO population, includi ng many Class I objects, in 



five nearby star forming regions. Similarly, 



Enoch et al.l fl2009[) identified 89 Class I YSOs 



using 1.25 fim to 1.1 mm data. I Whitney et al.l (120031 ) and iRobitaille et al.l (120061 ) showed 
how observables such as the SED, color, and polarization could be used to differentiate be- 
tween YSOs at diffe rent viewing angles as the y physically evolve. These results are the basis 
of a relatively new (Ivan Kempen et al.l 120091 ) classification sequence based on the physical 
characteristics of the YSO rather than its observational properties (e.g. SED, extinction, 
etc.). 

Although SEDs measured over a wide range of wavelengths are useful for general char- 
acterization of YSOs with their disks and envelopes, they do not probe in detail the region 
of the inner disk where the star and disk interact. Near-IR observations are sensitive to 
light from the protostellar photosphere, thermal emission from the warm (~ 1000 K) inner 
disk, as well as high excitation atomic and molecular transitions. Several of these transitions 
are used to probe the accretion onto the star (e.g. atomic hydrogen) as well as the wind 
and outflow (e.g. [Fe II] and H 2 ), whereas others are sensitive to t he temperature structure 
of the disk (e.g. CO, Calvet et al. 1991). iGreene &: Ladal (119961 ) presented the first large 
near-IR (1.15 /im - 2.42 /iin ) spectroscopic survey of pre-main sequence stars, including a 
number of Class I YSOs. They also observed FU Orionis-like stars, characterized by deep 
CO and water vapor absorption bands and without photospheric absorption lines (see section 
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7.4). They showed how extinction affects the appearance of the spectrum, and how gravity 
affects separate Class I YSOs and FU Orionis-like stars from Class II and HI objects by 
compa ring the equivalent widths of Na I + Ca I versus CO. The study by lReipurth fc Aspin 
( 119971 ) showed in all cases that the spectra of Herbig-Haro outflow sources are highly veiled, 
and in several cases the spectra are similar to FU Orionis-like stars. Recently, new instru- 
mentation has become ava ilable to obs erve the near-IR spectra over a wider bandwidth at 
higher spectral resolution. iBeckl (120071 ) used 5 different methods to estimate the extinction 
from the spectrum of 10 Class I and flat-spectrum protostars, resulting in widely different 
solutions, and she demonstrated that estimating the extinction for embedded protostars re- 
quires great care. She also estimated the spectral types of 5 protostars, and showed that 
the spectra of embedded protostars are charact erized by a strong water ice absorption band 
from 2.9 /im - 4.0 /im. iMuzerolle et al.l (119981 ) also noted the importance and difficulty of 
accurately estimating the extinction to highly extincted sources. 



Connelley et al.l (j2008al ) presented a new sample of Class I YSOs selected from the IRAS 



database to have increasing flux with wavelength, be nearby, and be optically obscured. We 
observed the near-IR spectra of Class I YSOs selected from this sample to address a number 
of questions. What are the spectral types of Class I YSOs? What fraction of Class Is are 
accreting, and h ow does this compare to T Tauri stars? What is the extinction and veiling 
of these objects? IConnelley et al.l ( j2008bl ) found that a large fraction of these protostars are 
binaries. Are the sp ectra of protostellar binaries related, as they often are for T Tauri stars 
(jMonin et al.ll2007| )? What can we learn about the interaction between the star and disk by 
observing the relationships between various emission lines? We discusses our sample selection 
and observations in §2. The atlas of Class I spectra is presented in §3. We discussion of 
specific spectral line features of circumstellar origin in §4. We discuss our method to estimate 
the spectral type of each protostar in §5, and in §6 we examine our methods of estimating 
extinction. We discuss our results in §7, and present commentary on selected targets in §8. 



2. Sample and Observations 



We observed a sample of Class I YSOs selected from the sample of IConnelley et al. 



(j2008al ). These objects were selected from the whole sky and not limited to previously 
known star forming regions. They were selected to have increasing flux density in Janskys 
with wavelength through the IRAS bands 0, to be in regions of nearby high optical extinction, 



1 Since spectral index is defined as -d(log^F^)/d(log^), it is possible for a source to have increasing flux 
density with wavelength yet still have a negative spectral index 
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and in most cases to have a red (H— K>1) near-IR counterpart observed by 2MASS. Since 
our observations were conducted at the NASA IRTF, we selected objects that are within the 
declination limits of the IRTF (5 < +70) and that rise above 2 airmasses f rom Mauna Kea 



(6 > — 40), which excluded 38 of 242 targets. In contrast to the study of IConnelley et al. 



( 12008af ). which is focused on imaging, spectroscopy requires brighter targets to get a high 
signal to noise observation in a reasonable amount of time. Our goal was to get a S/Nf^300 
spectrum for each target, which limited our selection of targets to sources brighter than 
K~ 12. The remaining sample included 124 targets brighter than K=12 and within the 
telescope's declination limits. We also observed targets that have an H-K color of ~ 1 or 
greater to ensure that each target is deeply embedded. Our K~ 12 magnitu de limit excluded 



many o f the deeply embedded, and thus very young, ob jects presented in IConnelley et al. 



(j2008al ) as none of the objects in IConnelley et al.l (j2008al ) with a spectral index greater than 
2.12 is brighter than K=12. 

We observed a total of 110 sources, which are listed in Table 1. This is the first time 
that infrared spectra of many of th ese sources have been pu blished. While these sources were 
all selected to be Class I YSOs by IConnelley et al.l (j2008al ). at least four are well known as 
T Tauri stars. Objects with a negative spectral index (Table 1) may be regarded as Class 
II (T Tauri) stars. The fact that there are a few well known T Tauri stars in this sample 
informs us as to the usefulness of these definitions near the ClassI/II boundary. Also, Class 
I state is observationally defined, and does not necesarily reflect the physical state of the 
YSO. Since the known T Tauri stars satisfied our selection criteria, they have been left in the 
analysis. When applicable, results are presented in Section 6 with these four T Tauri stars 
removed from the sample. Some of the objects presented here are the binary companions 
to Class I YSOs, and may not be Class I objects themselves. Higher angular resolution 
mid- and far-IR surveys will better distinguish embedded YSOs from nearb y T Tauri stars. 



Cautio n is advised as some targets have known close binary companions (IConnelley et al. 



2008a]) that are unresolved with the angular resolution of IRTF/SpeX, and light from both 
components are included in the spectum. Fifteen such objects are flagged in Table 1. An 
unknown number of targets have unresolved binary companions. 

The distribution of the spectral indices of these targets, based on IRAS data from 
12 /im to 100 /mi, is shown in Figure 1. The figure shows the spectral index distribu- 
tion for targets in the Taurus- Auriga and Perseus star forming regions in red, which has a 
mean spectral index of +0.39. The spectral index distribution for targets in the Orion and 
Ophiuchus star forming regions are also shown in gray, which has a mean spectral index 
of +0.64. These two distributions are more than 97.5% likely to be from different parent 
populations. The median spectral index of this sample of Class I objects is +0.50 a nd the 
average is +0.55. In comparison, the sample observed by IConnelley et al.l (j2008al ) had a 
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median spectral index of +0.79 and included targets with spectral indices greater than 3. 
The spectroscopically observed sample presented here has a lower median spectral index and 
includes no objects with a spectral index greater than +2.12. The bolome tric luminosities 



i n Tab le 1 were calculated using the IRAS fluxes using the expression in IConnelley et al. 



( 120071 ). In many cases there are multiple protostars within the IRAS beam, all contributing 
to the observed IRAS fluxes. As such, the calculated bolometric luminosity may be greater 
than the bolometric luminosity of the individual protostar that we observed. 

Each object w as observed with the 3.0 m IRTF on Mauna Kea, Hawaii with SpeX 



( IRayner et al.l 120031 ) in the short cross-dispersed mode, which covers 0.8 fim to 2.45 /im in 
each exposure. There is a gap in wavelength coverage between 1.82 /im and 1.88 fim, 
corresponding to a wavelength range where the atmosphere is relatively opaque. Several of 
the spectra in Appendix A are shown with a gap between 1.38 /mi and 1.41 fim where 
the spectra are typically very noisy due to the opacity of the atmosphere in that wavelength 
range. We used the 0'.'5 wide slit, which gives a resolution of R=1200. The star was nodded 
along the slit, with two exposures taken at each nod position. The individual exposure 
times were limited to two minutes to ensure that the telluric emission lines would cancel 
when consecutive images taken at alternate nod positions were differenced. An AO telluric 
standard star was observed after at least every other protostellar target for telluric correction, 
usually within 0.1 airmasses of the target. 

An argon lamp was observed for wavelength calibration and a quartz lamp for flat 
fielding. An arc/flat calibration set was observed for each target /standard p air. The SpeX 



data were flat fielded, extracted, and wavelength calibrated using Spextool (jCushing et al. 



2004 ) . After extraction and wavelength calibration, the individual extracted spectra were 
coadded with xcombspec. Xtellcor was then used to construct a telluric correction model 
using the observed AO standard, after which the observed spectrum of the target was divided 
by the telluric model. Finally, xmergeorders was used to combine the spectra in the separate 



orders into one continuous spectrum. These are all IDL routines written by ICushing et al. 



(120041 ) to completely reduce SpeX data. Spectral line flux, equivalent width, and FWHM 



were measured using the SPLOT routine in IRAF. 



3. Spectral Type Determination 

Among the most important properties of any star is the effective temperature or spectral 
type. We were able to measure enough photospheric absorption lines in 50 YSOs (45% of 
the sample) to be able to estimate the spectral type of the star based on these photospheric 
features alone. We used these results to estimate the photospheric spectral type distribution, 
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to determine if there is any trend of spectral type with spectral index, and if there is any 
regional dependence on spectral types. 

To estimate the photospheric spectral type of each protostar, we first measured the 
equivalent widths of 51 photospheric absorption lines (with uncertainties) in the spectrum of 
each Class I YSO. These lines were chosen so that some would have a high equivalent width 
for any given spectral type. We also measured the EWs of these 51 lines (wit h uncertainties 



as w ell) for several stars selected from the SpeX near-IR spectral library (j (Pushing et al. 



20051 ). Of these 51 lines, our program only considers lines in the YSO spectrum with a 
signal-to-noise ratio greater than a specified value, typically 2. To take into account the 
effect of veiling (excess continuum emission that appears to weaken spectral lines), we de- 
veiled the EW measurements from each protostar by using a grid of veiling temperatures and 
values of r^ (the ratio of the veiling flux divided by the continuum at K-band). We modeled 
the veiling emission as a single-temperature blackbody. To estimate the relative likelihood 
that the two EW measurements (one from the YSO and the other from the library star) 
match for a given line, we take the integral of the product of two Gaussian functions where 
the median of each Gaussian is the EW measurement and the standard deviation is the 
uncertainty in the EW measurement. This results in a relative "goodness-of-fit" parameter 
between each de-veiled line and each spectral standard line. We then compared the median 
value of this "goodness-of-fit" parameter for all of the lines for each of the spectral types 
we considered. Thus the spectral type with the highest median "goodness-of-fit" parameter 
is the best-fit spectral type for the protostar. If the goodness-of-fit parameter for a given 
spectral class is greater than half of the best-fit value, then we consider that spectral class 
within the reasonable range of possible spectral types for that YSO. The results of this 
spectral type matching are presented in Table 2. This procedure takes into account the S/N 
ratio of the observations, how precisely we could measure each line's EW, and the wavelength 
dependence of the veiling even though the amount of veiling or veiling temperature are not 
well known. Indeed, we use the result of this absorption line EW matching code to constrain 
Tk in the cases where we estimate the protostellar spectral type. We note that since the EW 
of a line is unaffected by extinction, the uncertainty in the extinction to each protostar does 
not directly affect the accuracy of our spectral type estimate. 

The results presented in Table 2 show that the uncertainties in the spectral type es- 
timates vary widely. This is caused by the spectra of the targets themselves rather than 
the code being unable to distinguish between stars of different spectral types. In order to 
test the code's precision, we gave the code the EW values for one of the Spex library stan- 
dard stars plus an error. This error was randomly chosen from a normal distribution with 
a FWHM of the EWs error bars as measured from the standard star's spectrum observed 
with SpeX. The code correctly identified the spectral type of the star with an uncertainty 



-7- 



of ±1 sub-class. A few of the spectral type estimates of Class I targets approach this pre- 
cision, but most do not. In some cases, the "goodness-of-fit" versus spectral type curve is 
double peaked. The precision of the match could be affected by flux from an unresolved 
binary companion. There are 15 known cases where light from both components of a close 
binary are included in the spectrum, but there are also an unknown number of cases where 
the target is an unresolved binary. Since our spectral type estimating routine only uses the 
EWs of narrow photospheric lines and does not use broad molecular features, circumstellar 
material could also affect the results of our matching only if that circumstellar spectrum 
has narrow "photospheric" lines characteristic of a different effective temperature than the 
central protostar. Also, some photospheric lines, such as from alkali metals, are known to be 



pressure sensitive ( IKirkpatrick et al.ll2006l ). The lines they observed in a young brown dwarf 
had EW values different than expected based on the effective temperature of the star. Thus, 
pressure sensitive lines may also reduce the precision of the spectral type matches for our 
Class I YSOs since protostars have lower gravity than the dwarf stars in the SpeX library. 

Although our spectral type estimates may have poor precision in so me cases, comparison 



with s pectral type estimates from other stu dies shows general agreem ent. IWhite fc Hillenbrand 



( 12004J ) used R=34,000 optical spectra and iDoppmann et al.l (120051 ) used R=18,000 spectra 



at 2 fim to estimate the spectral type of a number of embedded protostars. There are six 
targets where we can compare our spectral type estimate with the estimate from one of these 
studies. In all cases, the our spectral type estimate agrees with the estimate in one of the 
above mentioned papers within the mutual uncertainties. 

Figures 2 and 3 show two presentations of our sample's spectral type distribution. We 
created Figure 2 by simply averaging the goodness-of-fit versus spectral type results (nor- 
malized by the peak value for each star) for each YSOs where we were able to estimate the 
spectral type, creating the mostly likely spectral type distribution for the sample. Figure 
3 shows the histogram of the best-fit spectral types for each star. Both figures show that 
approximately half of the YSOs in our sample are M class stars, but there are a number 
of stars of earlier spectral types as well. The A class was a very poor fit for most stars, 
resulting in the low value in Figure 2, but this is the best-fit spectral class for 3 targets (all 
in Orion). Since we do not have near-IR B-class stellar spectral templates for comparison, 
we do not know if any of our targets would be best-fit by a B-class spectrum. Although no 
star had a best-fit spectral class from M5 to M7, most stars had a low but non-zero proba- 
bility of being M5 to M7, again resulting in a low value in Figure 2. The reason for having 
no best-fit spectral class later than M4 could be due to there being no proto-brown dwarfs 
in the sample, a possible effect of selecting our sample from IRAS (based on the Baraffe 
et al. (2002) models, a star with an M4 spectral type at 1 Myr age would have a mass of 
0.2 M Q and would have an apparent K magnitude without extinction of 10.1 at the distance 
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of the Orion star forming region, well within our K<12 apparent brightness limit). Another 
possible reason is that the photospheric line EWs for main sequence dwarfs in the M5 to 
M7 spectral classes may poorly match the EWs for protostars with this range of effective 
temperatures due to an effect such as lower gravity 

Class I protostars are expected to be descending the Hayashi track at a nearly constant 
effective temperature. As such, the effective temperature is largely set by the mass of the 
protostar. Half (21/42) of the stars with spectral type estimates have best fit spectral types 
from MO to M4 (T P jf from ~3800 K to ~3100 K). According to the models presented by 



Baraffe et al.l (119981 ) at an age of 1 Myr (likely older than our targets), these spectra l types 



correspond to masses from ~1.0 M to ~0.3 M Q . Models presented by lBaraffe et al.l (j2002f ) 
showed that the effective temperatures do not significantly change much in this mass range 
at ages earlier than 1 Myr, nor do they strongly depend on the initial gravity of the protostar. 

We observed that a number of targets have CO absorption in excess of what is expected 
from a M4 dwarf star, which is the latest class to be a best fit to any target in the sample. 
To determine if the excess CO absorption we observed was due to the star belonging to a 
high luminosity class (I, II, or III), we compared our observations with EW measurements 
of M-class dwarf and giant stars. In most cases, the EWs from the dwarf stars were the best 
fit. Figures 4 and 5 plot the equivalent width of Na I and Ca I versus the equivalent width 
of CO as a probe of the gravity of each star. Most of our targets lie near the dwarf locus on 
this plot, and away from the giant locus. The addition of veiling pushes the equivalent width 
measurements towards (0,0). Due to the location of our targets on this figure, we conclude 
that the excess of CO absorption is not due to the targets being from a high luminosity class. 
Section 6.4 discusses these objects further. 



4. Circumstellar Line Features 

This section disusses the observed properties of emission and absorption lines of circum- 
stellar origin. Several of these lines are association with the mass accr etion process, and thu s 



are correlated with Br 7 emission, a common mass accretion tracer (IMuzerolle et al.lll998l ). 
The relationships between Br 7 emission and other circumstellar lines are discussed in the 
sub-sections association with those lines. The frequency of common circumstellar emission 
and absorption lines are summarized in Table 3. The equivalent widths of the emission lines 
discussed in this section are tabulated in Table 4. 



The amount of veiling contributed by infrared excess emission from circumstellar ma- 
terial is seen to vary widely from object to object in this spectroscopic atlas. Many objects 
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show no evidence of veiling, whereas the veiling is so high in many others that no spectro- 
scopic features from the protostellar photosphere are apparent. Since the amount of veiling is 
difficult to precisely constrain, we have divided our sample into two veiling groups for much 
of our following analysis. Objects where the veiling is low enough to allow photospheric lines 
to be observable are described as "low veiling". Objects where the veiling is high enough 
so that no photospheric lines are apparent are described as "high veiling". The amount of 
veiling necessary to obscure photospheric lines is dependent on the photospheric spectral 
type since an A-type star has deep, broad hydrogen absorption lines that are easier to see 
despite veiling emission compared to the narrow metal lines in the spectrum of an M-star. 
For example, the highest veiling measured for an M-type star is r^ = 1.6^'g whereas the 
highest veiling measured for an early type star is r^ = 8.8I32. Since late-type stars are the 
most common in this sample, objects with high veiling are likely to have r^ > 5. Our proce- 
dure for determining spectral types and quantitatively measuring the veiling, when possible, 
are described in Section 3. 



4.1. Ca II Infrared Triplet 

Ca II has a well known near- infrared triplet line at 8498 A, 8542 A, and 8662 A. Ca II 



emiss ion is likely produced in the protostellar magnetospheric infall of gas ( lAzevedo et al. 



20061 ) and thus is a useful accretion tracer. If the lines are optically thick, the line flux 



ratios should be 1:1: 1 since these lines are very close in wavelength. Detailed modeling by 



Azevedo et al.l (120061 ) shows that the line ratios are typically 1.1-1.2, but can be as large as 
1.5. High resolution spectra of ~90 T Tauri stars, and found that the shortest wavelength 
line o f the triplet tends to have the highest peak and also tends to be the narrowest of the 



lines ( jHamann fc Perssonlll990l ). CitetWhi2004 presented high resolution spectra of this line 
for a large sample of T Tauri and Class I stars, showing the Ca II line profiles range from 
being very narrow to having a FWHM up to 250 kms -1 . 

Although only 20 Class I YSOs (18% ±4%) showed Ca II emission, 20/33 (61% ± 14%) 
targets with detectable flux in I-band have Ca II emission, making Ca II a very common 
(although rarely observed) emission feature for Class I YSOs. Being associated with mass 
accretion, and specifically magnetospheric infall, Ca II emission is strongly associated with 
HI Brackett (Br) 7. All targets that have Ca II emission also show emission from Br 7, 
but 20/25 targets with a detected continuum at I-band and Br 7 emission also have Ca II 
emission. Thus, ~20% of targets have Br 7 emission but no emission from Ca II. In terms of 
equivalent width (EW), Ca II is the strongest emission line observed from Class I protostars. 



We also found that the EW ratios of the Ca II infrared triplet can deviate significantly 
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from 1:1:1 line ratios. The mean ratio for EWs498 / EW 8 542 = 11, but ranges from 1.8 
to 0.8. Similarly the mean ratio for EW 8 498 / EW 8662 = 1.3, but ranges from 1.9 to 0.8. 



Gamiero et al.l (120061 ) observed a similar range of li ne ratios in the EWs o f the Ca II lines for 



DI Cep, showing that these ratios are also variable. lAzevedo et al.l ( 120061 ) note that variation 
in the gas temperature and mass accretion rate of their models do not completely account 
for this variability. 



4.2. He I 



The He I line at 1.0833 /mi is seen in both absorption and emissi on from many youn g 
stars. The line profile is sensitive to the kinematics of the stellar wind ( jDupree et al.lll992l ). 
Edwards et al.l (120061 ) presented observations of 38 He I line profiles from T Tauri stars 
observed at a spectral resolution of R=25,000, and demonstrated that young stars have a 
wide diversity of inner wind properties. In comparison, our data have much lower spectral 
resolution (R=l,200, £v~ 250 kms -1 ). Details of the line profiles of He I will be discussed 
in a future paper. 

The He I line was detected in 35 of our Class I YSOs. This comprises 32% of the whole 
sample, and 52% (35/67) of the targets where there was enough flux at J-band to detect 
the continuum. Among the He I lines detected, 69% (24/35) of the He I lines show sub- 
continuum absorption, and emission is seen in 77% (27/35) of them. Note that objects with 
He I absorption can also simultaneously have He I emission. We see deep He I absorption 
from the three FU Orionis-like stars where we detected the J-band continuum. Data from T 
Tauri stars show that the wind a ccelerating region is very close to the stellar photosphere, 
and is possibly the stellar corona (jEdwards et al.ll2006l ). However, the flux from FU Orionis- 
like stars is likely dominated by light from the disk (IHartmann fc Kenyonl Il985l ). yet we 
still see strong He I absorption, suggesting that in these cases the wind originates from the 
disk. The He I features for all of the FU Orionis-like stars in the sample show only strong 
blue-shifted absorption (i.e. no emission or non-blue-shifted absorption). Thus, we propose 
that this is a possible way to indicate if a candidate star is an FU Orionis-like star. However, 
some other stars also have a He I feature that shows only blue shifted absorption despite not 
sharing other spectroscopic features with FU Orionis-like stars, so this criterion cannot be 
used alone. 
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4.3. [Fe II] and H 2 



In this analysis we consider the [Fe II] line at 1.644 fim and the H2 line at 2.122 /mi. 
Both [Fe II] and H2 are well known tracers of winds via shock induced emission. Emission 
from these species is often seen far from the stellar source in Herbig-Haro flows. The excita- 
tion mechanism for H 2 within a few hundred AU of a young star remained poorly understood 
until recently. Possible excitation mechanisms include UV (either from the accretion flow 
onto the star or fr om the star itself), X-rays from the star's corona, or shocks from the wind. 
Beck et al.l (120081 ) found that the properties of the H2 emission (e.g. the emission morpholo- 
gies and spatial extent) are most consistent with shocked excitation from the wind rather 
than excitation of gas by radiation from the central star. [Fe II] emission is a very useful 
probe of win ds, especially in re gions where the extinction is high enough to preclude optical 
observations (jBally et al.ll2007l ). [Fe II] emission probes high excitation temp erature (11,000 



to 12 ,000 K) winds, and in particular fast (>30 kms x ) dissociative J-shocks ( IReipurth et al. 
2000h . 



The study by iBeck et al.l (120081 ) was based on observations of 6 Classical T Tauri stars. 
We seek to determine if H2 emission from Class I protostars is also shock excited. The [Fe II] 
and H 2 lines are also often observed in the spectra of Class I protostars, with 44/107 targets 
having [Fe II] emission and 47/110 having H 2 emission. Among the 44 Class I objects that 
show [Fe II] emission, 35 (80% ± 13%) also have H 2 emission. The probability that a random 
sampling of 44 spectra from this data set would yield 35 or more objects with H2 emission 
is less than 10~ 6 . Considering the strong correlation between the presence of [Fe II] and H2 
emission, we conclude that H 2 emission from Class I YSOs is also likely to be excited by 
shocks in winds. 

Figure 6 shows that there appears to be a weak trend of increasing [Fe II] equivalent 
width with increasing H2 EW. All targets with an H2 EW less (stronger emission) than 
-3 A also show [Fe II] in emission, provided there was enough flux to observe the H-band 
continuum. With two exceptions, all targets with [Fe II] EW less than -2 A have H 2 emission. 
With one exception, only targets with high veiling have an [Fe II] EW less than -7 A. 

Figure 7 shows that the Br 7 and H2 EWs appear to be related by a weak (the correlation 
coefficient is 0.03) inverse relationship. Targets with strong Br 7 tend not to have strong 
H 2 emission and vice versa. Only targets with high veiling have strong emission from H 2 
(EW<-7 A) or Br 7 (EW<-9 A). It may be possible that higher mass accretion rates onto 
the star suppress the H2 ex citation mechanis m. However, since H2 emission is believed to be 
shock excited in the wind ( IBeck et al.l 120081 ). which is accretion driven, this scenario seems 
unlikely. Another possibility is that the H 2 line flux is not dependent on the mass accretion 
rate, and that the higher veiling associated with higher Br 7 EWs reduces the observed H 2 
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EW. However, as stated above, only targets with high veiling have strong emission from H2, 
contrary to this scenario. Increasing veiling will tend to push the EW values towards (0,0) 
in Figure 7 (and in Figure 6 as well). Although many data points are near the origin of the 
figure, they are predominantly targets with low veiling. We also note that no early type star 
where Br 7 is seen in absorption shows H2 emission. Spectroscopic monitoring should show 
if the variability of the Br 7 and H2 lines are correlated or independent of each other. 

We observe that the line EWs for H 2 and [Fe II] are stronger in targets with high veiling. 
Objects with high veiling might be expected to have greater Br 7 luminosities since both 
are well known accretion tracers. Since the wind is powered by accretion, it would also be 
expected that targets with high mass accretion would have stronger winds. Howeve, if both 
Br 7 and veiling are accretion tracers, then why is there an inverse relationship between H 2 
and Br 7 EWs? These correlations show that the H 2 and [Fe II] lines, which are from the 
shocks in the wind, are not proportional to the instantaneous mass accretion rate (traced 
by Br 7). Since atomic hydrogen is only ionized very close to the star whereas H 2 traces 
colder more distant gas, the source of these lines are not co-located. Also, we speculate in 
section 6.4 that it may be possible that the mass accretion rate is quite high in the cases 
where there is strong H 2 and [Fe II] emission, but that the Br 7 emission mechanism has 
collapsed resulting in relatively lower Br 7 emission with increased mass accretion rate. We 
stress that this result is applicable to the region immediately around the central star, and 
not for an extended outflow. The width of our 0'.'5 slit corresponds to a range of ~70 AU at 
the distance of the Taurus star forming region and ~230 AU at the distance of the Orion 
star forming region. 



4.4. CO 

We have found that 23% ± 5% (25/110) of our targets sho w the CO band heads in 



emission. This is very close to the value for HH sources found by iReipurth &: Aspinl (119971 ) 
of 20%. Targets that show CO in emission tend to have high veiling such that no photospheric 
absorption lines are apparent. Of the targets that show CO in emission, 91% ± 20% (21/23) 
have high veiling whereas 9% ± 6% (2/23) have low veiling. Although 2 stars have CO 
emission and low veiling, these two stars are both early type stars, one of which has veiling 
that is low enough that the H I lines can be seen at shorter wavelengths but high enough 
to almost completely veil the Brackett series of lines (r^ = 8.8 for IRAS 05513—1024). 
Targets with CO emission more often have high veiling than the sample as a whole, for 
which 55% ± 7% (61/110) show high veiling. 

We have found that the CO and Br 7 features are closely related (the correlation co- 
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efficient of their EWs is 0.44). Figure 8 plots the EW of CO versus Br 7, and is divided 
into three parts. Region A includes all objects where CO is seen in emission. CO equivalent 
widths that are consistent with absorption from the photosphere of a dwarf star are included 
in region B, from early type stars with no CO absorption to M4, the latest type in our 
sample (see Section 3). The CO absorption in region C is in excess of what is expected to be 
observed from the photosphere of a dwarf star. This figure shows that targets with CO in 
emission also always show Br 7 emission (an important tracer of magnetospheric accretion) 
and almost always have high veiling (red symbols) as well. Figure 8 shows a general trend 
of greater CO emission with greater Br 7 emission, but with significant scatter. Br 7 is of- 
ten seen in emission when the CO absorption is consistent with a dwarf stellar photosphere 
(region B) and when the veiling is low (black symbols). There is a trend of increasing Br 7 
emission with decreasing CO band head absorption, possibly due to CO and/or veiling emis- 
sion increasing as the mass accretion rate increases. Notably, Br 7 never shows significant 
emission when the CO is in absorption beyond what is consistent with the star's spectral 
type. This situation is most commonly seen among FU Orionis-like stars (see section 6.4), 
which also have high veiling. 

We interpret this relationship between CO and Br 7 as follows. When CO and Br 7 
are in emission and the veiling is high, the accretion rate is quite high. The surface of the 



disk is hotter than the disk midplane, accounting for the CO emission (jNajita et al.lll996f ) 
and the lack of any absorption features from the disk. At lower accretion rates, Br 7 is still 
seen in emission. The veiling is low enough that photospheric absorption bands dominate 
the CO feature. The low veiling suggests that the warm (> 1000 K) emitting surface area 
in the disk is relatively low. In contrast, FU-Orionis like stars are believed to be experi- 
encing aburet_ofjnass accre tion. Strong CO absorption is observed from FU Orionis-like 
stars. iMuzerolle et al.l (119981 ) empirically derived the well-used relationship between Br 7 
emission and mass accretion rate with observations of T Tauri stars. Although FU Ori-like 
objects have ver y high accretion luminosities and thus very high mass accretion rates (up 



to 10" 4 Moyr" 1 flHartmannl 120001 ) ) . they lack many of the usual mass accretion signatures 
commonly observed from T Tauri and many other Class I YSOs (e.g. Br 7 or Paschen (3 emis- 
sion). Br 7 is not detected in some of the FU Ori-like objects (04073+3800, 18270-0153W, 
20568+5217, and 22051+5848) and only very marginally detected in the rest, as seen in 
region C of Figure 8. These observations show that the established relationship between 
observed Br 7 emission versus mass accretion rate is not valid for these objects. As such, 
FU Ori-like objects may accrete mass via a different mechanism than the magnetospheric 
funnel flows that produce the well known Br 7 emission in Class I and T Tauri stars with 
lower mass accretion rates. It is possible that Br 7 emission is proportional to the mass 
accretion rate at relatively low mass accretion rates, but the proportionality breaks down as 
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mass accretion rate increases to the point where at very high mass accretion rates there is no 
observed Br 7 emission. We postulate that the magnetosphere, which shapes the accretion 
flow for classical T Tauri stars, collapses when the accretion rate is as high as is found among 
FU Orionis objects. For further discussion of FU Orionis-like objects, see section 6.4. 

Figures 4 and 5 show the relationship between CO emission/absorption and Na I + Ca 
I. The FU Orionis-like objects in our sample lie in region E, below and to the right of the 
luminosity class III stars. A number of objects have sufficiently high mass accretion to push 
CO into emission, and these objects are located in region F of Figure 5. Absorption from 
Na I + Ca I is photospheric. However, Na I has also been observed in emission in many of 
these cases when veiling is high, resulting in occasionally negative values for the E W of (Na 
I + Ca I). Ca I was never seen in emission. 

CO emission and absorption trace gas at a moderate temperature (a few 1000 K) and 
high densit y (11 h > 10 9 c m" 3 ), however there are several potential excitation mechanisms. 
Figure 7 in ICalvet et all (119911 ) shows that young stars with a mass accretion rate below 
~ 10~ 7 MO yr _1 can have CO in emission due to stellar radiation heating of the disk surface 
or CO in absorption from the stellar photosphere. If this model is correct, then we should 
expect our sample to have many YSOs that are accreting mass (judging from their Br 7 
emission), have low veiling (suggesting a low mass accretion rate), with some YSOs showing 
CO in absorption and others showing CO in emission. Although we observed many examples 
of YSOs with low veiling and CO in absorption, CO emission is usually only observed when 
th e veiling is high and thus presumably when the mass accretion rate is high. Also, Figure 7 
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Calvet et al.l ( 119911 ) predicts that early type stars should have CO emission when the mass 
accretion is as high as 10~ 5 ' 5 M© yr _1 . Two of the eight early type stars in this study have 
CO in emission, whereas 6/8 have no CO emission or possible CO absorption (4 of these 
also have atomic hydrogen emission showing that they have a disk and are actively accreting 
mass). Both of the early type stars with CO emission also have high veiling, suggesting that 
the presence of CO emission is corre lated with high veiling and not with the spectral type of 
the central star. ICalvet et al.l ( )199ll ) also predicted that high mass ac c retion rates push the 
CO into absorption, as observed with FU Orionis-like stars. iMartinl (119971 ) proposed that 
the CO emission could be from the magnetospheric funnel flow of mass onto the star. If this 
hypothesis is correct, then multi-epoch observations should show that the CO band head 
flux should vary in step with other mas s accretion trace r s rela ted to the funnel flow, such as 
Br 7 or Ca II emission. More recently Idassgold et al.l ( 120041 ) concluded that CO emission 
arises in a thick layer near the surface of the disk where atomic hydrogen collisionally excites 
the CO. Velocity resolved observations can also help to determine the location of the source 
of the CO emission. 
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5. Extinction Estimation 

An accurate estimate of the extinction along the line of sight to each protostar is very 
important for deriving many properties of a YSO (e.g. luminosity, emission line fluxes, 
and properties derived from these). This is particularly challenging in the case of Class I 
YSOs that suffer from very high extinction, and where much of the observed flux may be 
scattere d light as shown by t he high fraction of candidate Class I YSOs with a reflection 



nebula (jConnelley et al .1120 071 ). Scattered light will cause a YSO to appear bluer, potentially 
leading to the extinction to be underestimated. Our estimate of the extinction starts with 
using a power l aw of the extinction ver sus wavelength based on empirical data. We used 



the results from iNishiyama et al.l ( 120091 ). who found that the J through K-band extinction 
is well fit by the power law A a oc A -2 . They also derived values for the extinction for the 
2MASS pass bands, allowing us to use 2MASS photometry to estimate extinction. Their 
result is based on observations of lines of sight towards the Galactic center. We note that 
the extinction law in star forming dark clouds may be different (e.g. Roman-Zuniga et al. 
2007). 

Extinction can be estimated in several ways. We used 2MASS JHK broad-band colors, 
modeling of the continuum, and emission line ratios to estimate the extinction to the targets 
in our sample. Not all methods could be applied to every target, and no method is perfect. 
Shocked emission lines from the outflow can be used, but they are seen though a different 
line of sight than the protostar even if the shocked emission is not spatially resolved. Our 
derived extinction values are summarized in Table 2. 



5.1. Broadband Colors 



We first used 2MASS J, H, and K band photometry of our targets to estimate the 
extinc tion to our targets. We used the values of Aa/A^- from Table 1 from INishiyama et al. 
(120091 ). Wi th these valu e s, we reddened the 2MASS colors of our targets to the T Tauri locus 
derived by iMeyer et al.l (119971 ). 2MASS did not detect some of our targets in J or H-band, 
so we were not able to estimate the extinction to all sources based on 2MASS photometry. 

Although this method of estimating the extinction is very simple, it has several caveats. 
2MASS observations cannot resolve close binary stars, so this method is not applicable to 
close binaries of comparable brightness. Scattered light from circumstellar material will 
make the object appear bluer and affect the derived extinction to the source. If an object 
is hidden behind an edge-on disk, then all of the observed flux may be scattered light. As 
noted above, it is very common for Class I YSOs to be associated with an infrared reflection 
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nebula. Table 1 flags objects where the K-band flux is likely to be dominated by scattered 
light (no point source is observed) or where scattered light could affect the near-IR colors 
(the object is associated wit h a reflection nebula) , a s judged from examining K-band images 
of their reflection nebulae in lConnelley et al.l ( 120071 ) . 



5.2. Continuum Modeling 

We used a simple procedure to model the observed spectrum of each protostar to esti- 
mate the extinction, veiling temperature, and the magnitude of the veiling. This program 
adds an infrared excess (i.e. veiling emission, which we simulate as a single temperature black 
body) to the spectrum of a star from the SpeX spectral library (Cushing et al. 2005, Rayner 
et al. 2008), then applies extinction to that sum using the reddening law described above. 
Naturally, the result of this procedure depends on the input spectral type. When possible, we 
used the best fit spectral type found by our equivalent width matching program; otherwise 
we used an M2V spectral type since that is the most common spectral type among the stars 
in this sample. Each model (star plus veiling behind extinction) is compared to the observed 
spectrum and the RMS of the fit is calculated. The code runs through a user-defined grid 
of veiling temperatures, veiling amplitude (r^), and visual extinctions (A„) to minimize the 
RMS difference between the observed spectrum and the model spectrum. We used veiling 
temperatures from 200 K to 2000 K, which is approximately the formation temperature of 
chondrules and refractory inclusions in meteorites (Boss & Grahm 1993, Alexander et al. 
2007). Contour plots of the RMS residual of the fit versus veiling temperature and extinc- 
tion show the confidence intervals for the fitted parameters. An example of the output of 
our modeling code is shown in Figure 9. Our 1 a uncertainties for veiling temperature and 
extinction are determined by the amount by which these parameters have to change from 
the best fit values for the RMS of the model to be twice the best fit RMS. In many cases, we 
removed strong emission lines (most often atomic and molecular hydrogen, and CO) from 
the observed spectrum to improve the RMS fit. 

The ability of our modeling code to precisely derive the extinction to the source varies 
greatly from target to target, and also depends on the spectral type of the template spectrum. 
This is reflected in the uncertainties presented in Table 2. We find that there is a degeneracy 
between the veiling temperature and the amount of extinction if the spectrum is dominated 
by a continuum that smoothly increases with wavelength, which is quite common among 
Class I YSOs. The amount of extinction is often very well constrained at a given veiling 
temperature, but good fits can often be found for a wide range of veiling temperatures and 
corresponding extinctions. We also find that among stars with late type spectra, the code 
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can often fit the observed spectrum using different amounts of veiling for different template 
spectra, and thus this procedure is a poor discriminator of spectral type. 

Figure 10 shows the histograms for extinction based on both near-IR broadband colors 
and continuum modeling. Figure 11 shows that the extinction estimates based on near-IR 
colors tend to be ~ 30% lower than the estimates based on continuum modeling. This 
discrepancy may be caused by the use of different photometric systems on very red objects, 
or scattered light from an associated reflection nebula affecting the observed colors more 
than the spectra. Another difference is that the 2MASS colors deredden to the T Tauri locus 
whereas the continuum modeling routine more accurately takes into account the contribution 
of the stellar photosphere and veiling emission without scattering. Although it would not 
account for this systematic offset, the large uncertainties in the extinction calculated from 
continuum modeling mean that these estimates are often consistent with the extinction 
estimates based on the near-IR photometry. The median value of the extinction distribution 
based on continuum modeling is 12.7 A^ magnitudes, whereas the median value of the 
extinction distribution based on 2MASS photometry is 9.8 A„ magnitudes. 

Our continuum modeling program also allowed us to derive a distribution for the veiling 
temperature, which peaks near 1500 K. This result is likely influenced by wavelength coverage 
of our data, since the peak of a 1500 K blackbody is at 1.9 ^m, near the middle of our 
wavelength coverage. 

We found it very difficult to precisely constrain the amount of veiling (r fe ) for most of 
the stars in the sample. For objects whose spectrum is dominated by veiling, we can only 
say that the veiling is high, likely greater than r^ ~ 5 since r&=8.8 is the highest veiling value 
that we were able to meaningfully constrain. In the cases where the spectrum shows clear 
photospheric lines, these lines are used to constrain the veiling. In these cases, uncertainty in 
the estimate of the spectral type of the star limits the uncertainty of the veiling. Therefore, 
for much of the analysis regarding veiling, we have simply divided this sample into two 
groups: targets with "high" veiling and targets with "low" veiling, as defined in section 4.0. 

5.3. Emission Line Ratios 

If the intrinsic line ratio of a given pair of emission lines is known, the observed line 
flux ratio can be used to calculate the extinction to the source. Lines that have the same 
upper state are particularly useful since the line flux ratio only depends on the transition 
probabilities to the two lower states. In this case, there are no assumptions about the 
emitting gas being in LTE or concerns that the two lines may be being emitted by gas with 
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different properties (temperature, density, distance from the source, etc.). In other cases it 
is reasonable to assume LTE, and thus the intrinsic line ratio can be calculated. 



[Fe II] emission lines at 1.644 /im and 1.257 /im are often seen in the spectra of 
young mass accreting stars. The ratio of these line flu xes has previously been used to 
estimate the extinction to the Cas A s upernova remnant (lEriksen et al.ll2009l ) and to make 
an extinction map towards an AGN (jStorchi-Bergmann et all 120091 ) since these two lines 
also share the same upper level. iReipurth et all (120001 ) also discussed the use of these lines 
to estimate extinction. Since these lines are forbidden, they should be optically thin. As 
such, this pair of lines is a useful tool to estimate the extinction to a win d source. Th e 



i ntrins ic line ratio is not yet well constrained, as noted by lEriksen et all (120091 ) . lEriksen et al. 
(120091) used a value o f [Fe I I] t o R7 / [Fe H] 1.^44. = 1.49, based on observations of P Cygni 
bv ISmith fe Hartiganl (hood). IStorchi-Bergmann et all (120091 ) used a value of 1.36 based 
on ca lculations (INussbaumer fe Storeylll988l ) verified by observations (IBautista fe Pradham 
19981 ). For the purpose of this paper, we adopt a value of [Fe II] 1.257 / [Fe II] 1. 644 = 1.36. 



We observed both the 1.257 /im and the 1.644 /im [Fe II] lines in 25 targets. Our 
extinction estimation used the same reddening law as with all of the other methods we used. 
Generally, our extinction estimates based on the [Fe II] line ratios are consistent with the es- 
timates based on continuum modeling, taking into account the mutual uncertainties. In the 
case of IRAS 04286+1801, the [Fe II] lines were the only reliable way to estimate the extinc- 
tion to this target among the methods we used. In this case, the near-IR flux is dominated 
by scattered light, making the object appear bluer and thus leading to an underestimated 
extinction. Furthermore, being an FU Orionis-like object, there is no apparent flux from 
the photosphere or atomic hydrogen emission (see Section 6.4). Since the [Fe II] emission 
is from the wind, the extinction calculated from the [Fe II] line ratio may be different from 
the true extinction to the protostar itself. However, considering that the slit was 0'.'5 wide, 
the [Fe II] emission in our spectrum is very near the location of the protostar, and not from 
a distant part of the outflow. We note that the 1.644 /im was occasionally on the edge of a 
comparably strong Br 12 absorption line or emission line. Such cases are designated with an 
asterisk in Table 2. Higher spectral resolution would be useful in separating these features 
to get a better flux measurement of the [Fe II] line. 

The Br 7 (2.16612 /im) to Paschen f3 (1.28216 /im) line ratio can be used as an extinction 
tracer assuming Case B recombination. In this case, the intrinsic line ratio is 5.75±0.15 
(jStorey fe Hummer! I1995T ). Although most of our targets show some Br 7 emission, the Pa 
13 emission line was detected in only about a third of them, usually because of high J-band 
extinction and as such there being no detectable continuum at J-band. We did not extract 
orders for which there was no continuum flux, so in these cases we did not measure the Pa (3 
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equivalent width. In those cases where both lines were measured, we used them to calculate 
the extinction to the target, and these values are listed in Table 2. 

In the case of em bedded prot ostars, the H 2 v=l-0 Q(3) (2.42373 /mi) to S(l) (2.12183 fj,m) 



lines have been used (IBeckl 120071 ) since these two lines share the same upper state. Many of 
our Class I YSOs have strong H 2 emission, and we attempted to use these lines to estimate 
the extinction to these sources. We found that the extinction derived from these H2 lines 
often greatly differed from the extinction estimates based on broad band colors or continuum 
modeling, and was occasionally negative. We determined that this line ratio is an unreliable 
extinction estimator since in 4/10 cases, this line ratio was observed to be variable. By 
coincidence, the v=l-0 Q(3) line is very close to a narrow (5\ = 7.1 A) but opaque telluric 
absorption line at 2.42412 /im (in vacuum). The Doppler shift due to the Earth's orbital 
motion moves this telluric line such that the v=l-0 Q(3) line can be either on the blue wing of 
the line (and thus the Q(3) line would be partially absorbed) or completely clear of it. This 
telluric feature as well as the H2 lines are not resolved at the spectral resolution (R=1200) 
of our data. Thus, a significant amount of the flux of the v=l-0 Q(3) line may be lost to 
this telluric absorption feature and not recovered when the telluric correction is carried out 
during the data reduction q This problem may be ameliorated by using sufficiently high 
spectral resolution to resolve this telluric absorption feature to accurately compensate for 
it. In the cases where the v=l-0 Q(3) to S(l) line ratio did not vary, the observed line 
ratio was used to estimate the extinction with the understanding that the extinction may 
be underestimated if a significant amount of the flux of the Q(3) line is lost to this telluric 
absorption feature. 



6. Discussion 



We now discuss several results revealed in the preceding analysis. We particularly 
examine the physical implications for our observed sample of embedded protostars. 



2 Consider a telluric absorption line that is narrow (and thus unresolved) but very deep. If the telluric 
line transmits 10% of the incident flux, then the telluric correction at that wavelength should be a factor of 
10. However, since the line is unresolved, the observed depth of this line will be much less (for example, 70% 
transmission) and the telluric correction will be similarly less (1.4 rather than a factor of 10). 
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6.1. Effect of Nebulosity 

Having noticed that many of the stars with FU Orionis-like spectra (see section 6.4) were 
associated with reflection nebulae, we investigated the connection between reflection nebu- 



l osity and the near-IR spectrum. We used the K-band images presented in IConnelley et al. 



( 120071 ) to classify the nebulae into three groups: strong nebulosity (no K-band point source 
is seen), some nebulosity (a K-band point source is seen along with a reflection nebula), or 
no nebulosity. The reflection nebulae have an average size of ~ 10000 AU or ~20" on the 
sky. The spectra as a whole were divided into three parts: spectra with "high" veiling (no 
photospheric lines apparent), spectra with "low" veiling (the spectrum shows photospheric 
absorption lines), and spectra that are similar to FU Orionis-like stars with deep water and 
CO absorption and without photospheric absorption lines. 

The results of our analysis are gathered in Table 5. We found that although there are 
only 10 targets with FU Ori-like spectra (9% of the sample), 40% of the YSOs with strong 
nebulosity have this type of spectrum. Furthermore, all objects in our sample with FU Ori- 
like spectra are associated with a reflection nebula. An exception is IRAS 06297+1021W, 
which we consider to be an FU Orionis-like object despite having emission lines uncharac- 
teristic of such an objct, and which is not associated with a reflecton nebula. YSOs with 
strong nebulosity have a higher fraction of FU Ori-like spectra (6/15) than objects with- 
out nebulosity (0/54) with greater than a 99.9% confidence. We found that the fraction of 
objects with high veiling has no dependence on the presence or absence of a reflection neb- 
ula. However, objects with low veiling are less likely to be found among targets with strong 
nebulosity than objects without nebulosity with a 99.9% confide nce. These binomial co nfi- 



dence intervals were calculated using the expressions presented by lBrandeker et al.l ( 120061 ) in 
their appendix B2. Targets without nebulosity and those with some nebulosity have similar 
fractions of objects with high veiling, low veiling, and FU Ori-like spectra, so these two 
nebulosity groups are difficult to distinguish based on their spectra. In summary, targets in 
this sample with strong nebulosity are less likely to have low veiling and more likely to be 
an FU Orionis-like star than targets with some nebulosity or no nebulosity. 

These relationships show that there is a correlation between the appearance of the 
YSOs on large scales (> 1000 AU) and the properties of the star and inner disk at very 
small scales. Objects with strong nebulosity could be objects with edge-on disks that block 
the central star from view, in which case we would see no K-band point source and all of 
the light we see is scattere d off of circumstellar materia l. This is indeed the case for many 



of the objects observed by IWhite fc Hillenbrand! (120041 ). who observed objects apparently 



more embedded than T Tauri stars (many objects were only seen in scattered light) and 
observed photospheric features in all of them. If objects with strong nebulosity are merely 
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seen edge-on, then we would expect to observe the same fraction of objects with low veiling 
or FU Orionis-like spectra in all of the nebulosity groups. However, this is not what we 
observe. In order to not see a K-band point source, these objects must therefore be more 
deeply embedded on average, and presumably less evolved, than objects with some or no 
nebulosity. Among the objects with no reflection nebula, there may be no objects with FU 
Orionis-like spectra because either objects with no nebulosity are too evolved or because the 
increase in the luminosity of the object in the FU Orionis phase would illuminate nearby 
circumstellar material, creating a reflection nebula that previously was not visible. 

Why are there so few objects with low veiling among objects with strong nebulosity, 
and so many objects with low veiling that have some or no nebulosity? The circumstellar 
material responsible for veiling also scatters and absorbs light, and can also obscure the 
central star. Objects with low veiling are most common among targets without a reflection 
nebula and quite rare among targets with strong nebulosity. If objects with low veiling were 
simply less luminous and thus do not illuminate a reflection nebula as often, then we would 
not see an equal number of low veiling objects among targets with some nebulosity and 
without nebulosity. One possible explanation is that objects with some or no nebulosity are 
simply more evolved, with less circumstellar material in the envelope (accounting for the 
optically thinner cloud, allowing us to see the central star) and less material in the inner 
disk (accounting for the lack of veiling emission). This would then suggest that FU Orionis 
objects, which are most common among targets with strong nebulosity, would be relatively 
less evolved. 



6.2. Spectral Type Dependence on Star Forming Region 



The spectral type distributions of pre-main sequence stars has be en shown to have a 



strong dependence on the star forming region (ILuhman et al.ll2003l ). iHillenbrandl (119971 ) 
compared the spectral type distributions for the Ophiuchus, Chamaeleon, Taurus/ Auriga, 
Lupus, L1641, and Orion Nebula Cluster (ONC) star forming regions. Both the ONC and 
the Ophiuchus star forming regions have a number of early type stars, extending to OB 
in Orion and B type in Ophiuchus. However, the spectral type distribution of the ONC 
peaks at M3 whereas the spectral type distribution of Ophiuchus peaks earlier than M-class. 
The spectral type distribution for the Taurus star forming region (e.g. Briceno et al. 2002, 
Luhman et al. 2003, Luhman et al. 2006) peaks near K7, with most stars having spectral 
types from K0 to M8. While these studies found no stars i n the Taurus star forming region 
with a spectral type earlier than G, iKenyon fc Hartmannl ( 119951 ) found 3 young stars in the 
Taurus/Auriga clouds with spectral types earlier than G. 



-22 - 



To determine if our sample of Class I YSOs shows a spectral type dependence versus 
star forming region similar to the pre-main sequence studies described above, we divided 
our sample into two groups of roughly equal size. Group 1 consists of targets in the Taurus, 
Auriga, and Perseus star forming regions, and has 21 objects. These clouds are nearby, 
have a low stellar density and are not forming massive stars. Group 2 consists of targets in 
the Orion and Ophiuchus clouds, and also has 21 objects. The distributions of the best-fit 
spectral types are shown in Figure 3. The spectral type distributions are overall very similar. 
The two-sample K-S test shows that the best-fit spectral type distributions have less than a 
90% chance of being drawn from different parent populations. However, Group 1 has no star 
with a best fit spectral type earler than F5. In contrast, 6/21 (29%) of the stars in Group 2 
are A or F stars. If the 4 known T Tauri stars are excluded, it remains true that Group 1 
still has no star with a spectral type earlier than F5. Group 2 would then have 5/20 (25%) of 
stars with an A or F spectral type. The regions in which we found A and F type stars (Orion 
and Ophiuchus) are associated with well known populations of early type stars, in this case 
the Orion Nebula cluster and the Sco-Cen OB association. We stress that these results are 
based on those stars for which we could determine a spectral type, and are only applicable 
to the clouds as a whole assuming that protostars of each spectral type are equally likely 
to have low enough veiling to allow a spectral type to be estimated. We also note that the 
Serpens clouds have very few targets for which a spectral type could be estimated. Of the 
eleven targets in that region, the spectra of only two targets have photospheric absorption 
lines, and we could get a spectral type estimate for only one target. 



6.3. Class I Accretion Fraction 



An important question to address is the fraction of Class I YSOs that are actively 
accreting mass. Unfortunately, the fraction of mass accreting T Tauri stars has not been 
well characterized, in part because the result depends on the age of the T association. The 
fraction of T Tauri stars with IR excesses, and thus warm disks, drops from ~50% to ~3% in 
a survey of sta r clusters ranging in age from 2.5 to 30 Myr (Haisch et al. 2001a, Hernandez 
et al. 2008). lHaisch et al.l (12001bh found that 65% of young stars in IC 348 have disks, 
whereas the younger NGC 2264 and Trapezium clusters have disk fractions of 86% and 80%. 
Luhman et al.l (120101 ) found that the disk fraction of the Taurus star forming region (N(Class 
II)/N(Class II and III)) i s ~75% based on 34 8 objects observed by Spitzer. Using this metric, 
the results presented by lEvans et al.l (120091 ) show that the disk fraction in Chamaeleon II 
is 78%±18%, in Lupus is 66%±9%, in Perseus is 89%±6%, in Serpens is 81%±7%, and in 
Ophiuchus is 81%±6%. The presence of a disk does not mean that all of these stars are 
accreting matter, but only that they could. 
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We used the presence of Br 7 emission (Muzerolle et al. 1998) or high veiling (when the 
veiling is high enough that no photospheric absorption lines are observed, r^ > 5, character- 
istic of a warm inner disk close to the star) to determine the fraction of our targets that are 
accreting matter. 93 stars show Br 7 emission, whereas 8 more have high veiling but no Br 
7 emission (these are mostly FU Orionis-like objects). Thus 92% ± 9% (101/110) of Class I 
YSOs are actively accreting matter at any given time. Excluding the 4 known T Tauri stars, 
the fraction of Class I YSOs accreting matter remains 92% ± 9% (98/106). 

The histogram of Br 7 equivalent widths (EW) in Figure 12 has been divided between 
targets with high and low veiling. Targets with high veiling have a mean EW of —4.21 
and a median EW of —2.93, whereas targets with low veiling have a mean EW of —1.20 
and a median EW of —1.63. Equivalent widths are measured relative to the continuum 
flux, which is elevated by the veiling in the case of targets with high veiling. Thus, targets 
with high veiling on average have much higher Br 7 line fluxes and presumably high mass 
accretion rates as well. For c omparison, the Br 7 E W histogram for the T Tauri stars 



with Br 7 emission observed by iMuzerolle et al.l (119981 ) has been overlaid in Figure 12. The 



two-sample K-S test shows that there is less than a 90% chance that the T Tauri Br 7 
EW distribution is drawn from a different parent population than Br 7 equivalent width 
distribution of the whole sample of Class I YSOs. Thus, the Br 7 EW distributions for 
accreting T Tauri stars and accreting Cl ass I protostars are indistinguishable with these 



data. We note that IMuzerolle et al.l (119981 ) selected T Tauri stars with previously measured 
mass accretion rates, so the targets they selected are not representative of T Tauri stars as 
a whole considering the large number of non-accreting weak line T Tauri stars. Although 
the Br 7 E W distributions are cons istent, Class I YSOs generally have higher veiling than T 



Tauri stars (IDoppmann et al.ll2005l ). so therefore Class I YSOs likely have a higher average 



Br 7 line flux and presumably higher mass accretion rates. 

There are 9 targets in this sample without Br 7 emission or high veiling, and thus do 
not appear to be currently accreting matter. We have estimated the spectral types of 8 of 
these. Among the non-accreting Class Is, early type stars are over represented and late type 
stars are under represented. The two-sample K-S test shows that there is a 90% chance that 
the spectral type distributions for accreting and non-accreting Class I YSOs are drawn from 
different parent populations. Three of these 9 stars are A or F type stars. An infrared excess 
or Br 7 emission that would be readily apparent for a late type star may be overwhelmed 
by the greater photospheric flux from these stars. Thus, these apparently non-accreting 
stars may still be accreting some matter. All targets, including the apparent non-accretors, 
were selected to have increasing flux with wavelength as observed by IRAS, so they all have 
significant circumstellar material. Non-accretors (excluding targets where there was another 
YSO in the IRAS beam) also tend to have a higher median spectral index than the sample 
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as a whole (1.3 vs. 0.5 from A =12 fim to 100 /xm), suggesting that the non-accretors tend 
to have relatively less warm dust and more cool dust than the accretors, consistent with the 
scenario of a cleared disk hole. 

We expected that the accretion luminosity (Lboiometric - L star ) would be proportional to 
the veiling. Accretion luminosity could only be estimated for those stars for which we have 
estimated the photospheric spectral type. Analysis of the bolometric luminosity (determined 
from IRAS observations) or the accretion luminosity showed no trend with our derived values 
of the veiling (r^). There are several objects with a late type spectral classification, very high 
accretion luminosities, and negligible veiling. For many of these cases there are other deeply 
embedded YSOs within the broad IRAS beam. To show how the veiling is related to the 
accretion luminosity, high angular resolution mid- and far-IR data are needed to spatially 
resolve these YSOs and higher spectral resolution near-IR spectra may be needed to more 
accurately determine the photospheric spectral types. 



6.4. FU Orionis-like Stars and Excess CO Absorption 

Several targets have spectra that cannot be well modeled by adding extinction and veil- 
ing to a stellar photosphere. The most common of these are targets with FU Orionis-like 
spectra, characterized by having a near-IR spectrum dominated by water vapor absorption, 
without clearly defined narrow photospheric absorption lines, few if any emission lines, and 
deep CO absorption in excess of what is observed in the spectra of late type dwarf pho- 
tospheres. Flux from these objects is dominated by emission from the disk (presumably 
due to accretion luminosity) with negligible flux from the stellar photosphere (Hartmann 
& Kenyon 1996). Since the heating of the disk is dominated by viscous dissipation versus 
stellar irradiation, the disk interior is hotter than the surface. Cooler material above the 
disk photosphere imprints the observed water and CO absorption bands on the spectrum. 
One object (04073+3800) has a highly veiled s pectrum without excess water absorption and 



has been identified as a FU Orionis-like star ( ISandell &: Aspinl Il998l ). FU Orionis-like ob- 
jects are found in region E of Figures 4 and 5. As a group, these targets are much more 
luminous than the rest of the observed sample. The median bolometric luminosity for tar- 
gets with FU Orionis-like spectra is 28 L whereas the median bolometric luminosity for 
targets in region D (with dwarf like spectra and gravities) of Figures 4 and 5 is 4.0 L Q . The 
two-sample Kolmogorov-Smirnov test shows that those targets with FU Orionis-like spectra 
have a higher bolometric luminosity than those targets with dwarf like spectra with 99.5% 
confidence. Although the bolometric luminosities are higher, the spectral index distributions 
(both of which are measured using IRAS fluxes from 12 fim to 100 /im) for region D and E 
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are statistically indistinguishable. Since these targets have deep CO absorption and tend to 
be excessively luminous, these are likely to be examples of FU Orionis-like objects. Indeed, 
two targets in our sample with this kind of spectrum, IRAS 04287+1801 (L1551 IRS5) and 
IRAS 21454+4718 (V1735 Cygni) are well known FU Orionis-like objects. 

Of the ten FU Ori and FU Ori-like stars in our sample, two (and possibly a third) 
are newly identified FU Ori-like candidates and these are the first near-IR spectra of them. 
FU Orionis-like objects are relatively rare. Previously, ~21 have been identified (Reipurth 
& Aspin 2010) and only 9.1% ± 2.9% (10/110) of the sample targets have this type of 
spectrum. Thus, the fraction of Class I YSOs that pass through a FU Orionis phase times 
the average duty cycle of the FU Orionis phase is roughly 9%. It is significant to find such 
a large fraction of FU Orionis-like stars among a sample of Class I YSO. Identifying two 
new embedded FU Orionis-like stars is particularly important since these observations have 
significantly increased the number of known embedded FU Orionis-like objects. Previously, 
only ~ 8 have been identified. Table 6 lists the FU Ori and FU Ori-like stars in our sample, 
as well as other stars that show excess CO absorption. 

In addition to FU Orionis-like stars, seven other stars show CO absorption greater 
than what would be expected from the star's spectral type. All of these show photospheric 
absorption lines, and in four cases the CO EW is only marginally greater than the M4 value 
and we do not consider these four to be candidate FU Orionis-like objects. IRAS 06393+0913 
has strong CO absorption but lacks the high veiling and deep water absorption bands that 
are commonly observed in the specrta of FU Orionis-like objects, and thus we do not classify 
it as an FU Orionis-like object. Three targets (IRAS 05327-0457W, 05404-0948, and 
16235—2416) have an A or F-type photosphere with weak CO band-head absorption. IRAS 
05404—0948 also appears to have weak water absorption bands. In these cases, the CO band- 
head and water absorption may be caused by circumstellar material or an unresolved late 
type companion. The absolute K-band magnitudes (simply the observed K-band magnitude 
of these objects corrected for distance, and not corrected for accretion or extinction) of the 
K and M-type stars in this sample are on average 2.5 magnitudes fainter than the absolute 
K-band magnitudes of the A and F-type stars (in comparison, the bolometric luminosities 
of main-sequence AO and M0 stars are different by ~ 8 magnitudes). However, the brightest 
K and M-type stars have the same absolute K-band magnitude as the faintest A and F-type 
stars, most likely due to a combination of accretion and extinction related effects. As such, 
for Class I YSOs it is reasonable to expect to see a contribution from a late type companion 
star in the spectrum of an early-type primary star. The other three stars (04292+2422W, 
04591-0856, 18341-0113N) are best fit with G or K type spectra but have CO absorption 
in excess of the latest type star that is a good fit. Remarkably, these six objects also tend to 
not have emission lines. IRAS 04591—0856 has H 2 emission lines, which may be associated 
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with a wind (IBeck et al.ll2008l ). and 05327— 045 7W shows weak Br 7 emission at the bottom 
of that deep absorption feature. Considering that these six objects share some properties in 
common with FU Ori-like stars (excess CO absorption, a dearth of emission features), these 
objects may be experiencing weak FU Orionis-like activity. Thus the shape of the continuum 
and the depth of the CO band-heads in the spectrum of a Class I protostar may not always 
be representative of the true spectral type of the stellar photosphere. 



6.5. Herbig Ae/Be Stars 



Herbig Ae/Be stars (lHerbigill960l ) were identified as stars with spectral type earli e r than 
F0 with Balmer emission lines and associated with a dark cloud. iHillenbrand et al.l (119921 ) 
showed that these are intermediate mass pre-main sequence stars with massive accretion 
disks. This spectroscopic survey of Class I YSOs found 3 targets with an A class best-fit 
spectral type (see Table 2), as determined by our spectral line matching routine. Since these 
objects satisfied our criteria to be Class I YSOs, and since the spectra of two of these stars 
show emission features, all three show veiling, we believe that these are embedded Herbig 
Ae stars. In this case, only 2.7% ± 1.6% of the Class I YSOs are embedded Ae type stars. 



6.6. Spectra of Binaries 



Many T Taur i binary systems have been spectroscopically observed, and Table 1 in 
Monin et al.l (120071 ) lists the number of T Tauri binary systems with classical T Tauri (pos- 
sessing a disk and accreting matter) and weak line T Tauri (non-accreting) components 



We used the number of T Tauri binaries in the center column of Table 1 in Monin et al. 



( 120071 ) to exclude objects with very clos e comp anions, and excluded the three systems with 
a passive disk. Table 1 in iMonin et al.l ( 120071 ) includes 59 (74.6% ± 9.7%) T Tauri binary 
systems where both components are classical T Tauri stars or both are weak line T Tauri 
stars, whereas 20 (25.3% ± 5.7%) are mixed systems. Has this correlation been established 
by the Class I phase? 

Over the course of this survey, we observed the spatially resolved sp ectra of both of 



the com ponents of 9 previously identified multiple protostars included in IConnelley et al. 



( l2008a[ ). This includes 1 triple system, for a total of 10 pairs. In order to be clearly resolved, 
these objects are well separated on the sky, with a median angular separation of 4'.'5. These 
objects have a median projected separation of 2700 AU, with a minimum separation of 
400 AU and maximum separation of 4500 AU. Both components of each pair also had to 
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be brighter than K=12 to get a high S/N observation. The properties of these binaries are 
given in Table 7. 

We noticed that the spectra of Class I binaries tend to have the same veiling, i.e. 
both components have high veiling or both have low veiling. Are the veiling of the binary 
components consistent with randomly pairing spectra from the whole sample? In the whole 
sample, 57.4% of targets have high veiling, whereas 42.6% have low veiling. Among the binary 
spectra, both components of 8/10 pairs (80% ± 28%) either both have high veiling or both 
have low veiling. The probability that less than 3 out of 10 pairs would have different veilings, 
having been randomly selected from this atlas of protostellar spectra, is 4.6%. Thus, we can 
state at the ~95% confidence level that binary protostars are not randomly paired and that 
their spectra tend to have the same veiling. Excluding the FS Tau system, which includes a 
well known T Tauri star, both components of 8/9 pairs (89% ± 31%) either both have high 
veiling or both have low veiling. Subjectively, the spectra of binary components often look 
remarkably similar. For example, the spectra of the components of IRAS 18383+0059 and 
the components of IRAS 05375—0040 are quite difficult to tell apart. We have shown that 
the trend of T Tauri binary components to both be accreting (or not) extends to the Class I 
phase at much higher accretion rates, for targets selected across the sky. We note that 8/11 
of the components with low veiling show Br 7 emission, and thus are accreting. 

Interactions between the circumstellar disks could explain the dearth of mixed Class I 
binary systems. For example, a close passage between two stars could trigger higher accretion 
and/or higher veiling in both systems. The projected separations between the stars of these 
binaries is within an order of magnitude of the expected size of a circumstellar disk, so 
interactions between the disks are plausible. If this scenario is true, then we would expect 
that closer binaries would tend to have higher veiling. Based on the values in Table 7, there 
is no trend of increasing veiling with decreasing projected separation. 

If interactions between the disks are not responsible for these stars having the same 
veiling, and if we assume that these protostellar binaries are coeval, then another explanation 
is that the veiling systematically changes with time within the Class I phase. Presumably, 
younger objects have higher veiling and more evolved objects have lower veiling. If the 
veiling of a given protostar were not to systematically change with time, but rather changed 
randomly on a timescale much less than the Class I life time, then the number of binary 
pairs with the same (or different) veilings would be indistinguishable from a random selection 
from the whole sample. 



H 2 emission is seen from 42.7% of the targets in the whole sample. Similarly, iDoppmann et al 

( 120051 ) found that 44% of Class I and flat spectrum YSOs have H2 emission in their smaller 
sample. Among the binary spectra, both components of 9/10 pairs (90% ±30%) either have 
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H2 emission or both do not have H2 emission. The probability that only 1 out of 10 pairs 
would have different H 2 emission, having been randomly selected from this atlas of protostel- 
lar spectra, is 2.8%. The presence of H 2 emission in both components of protostellar binaries 
is correlated with ~97% confidence. We also note that only 5 of 19 binary components 
(26.3%) have H2 emission, less than the 42.7% for the whole sample. 



6.7. Gravity and Late Type Continuum Profiles 



Lucas et al.l (120011 ) found that a number of young (~1 Myr) low mass (< 0.1 M ) 
sources in the Trapezium cluster have an H-band continuum spectrum that is characterized 
by a triangular profile with a peak near 1.68 /mi. They attribute this to the expected lower 
gravity of younger sources, and that this is an independent way to verify that these are 
substellar sources. A field dwarf observed by iKirkpatrick et al.l ( 120061 ) also has a triangular 
shaped H-band continuum spectr um which they inter pret as a sign of low gravity and thus 
youth since regular field dwarfs ( jCushing et al.l 120051 ) do not exhibit this feature. Several 
other young brown dwarfs have been identified (e.g. Gizis et al. 2002, Looper et al. 2007) 
with this feature. The surface gravities of protostars an d young brown dwarfs are very 
similar. The average log(g) for protostars from Table 3 of Doppmann et all ( 2005}) is 3.98, 
whereas average log(g) for 5 Myr brown dwarfs from Table 4 of iRice et al.l ( 120101 ) is 3.87. 
The triangular shape of the H-band continuum spectrum is caus ed by the reduction of the 
pressure induced opacity in H 2 in the atmosphere of the star ( iMarley et al.l l2010h . The 
spectrum of pressure induced opacity in H 2 is very broad, which has the effect of smoothing 
the opacity spectrum of the stellar atmosphere as a whole. Without this opacity source, the 
shape of the continuum near H-band is dominated by the opacity spectrum of water. 

We used our continuum fitting routine described above to determine which objects in 
our sample with a late type continuum (excluding FU Orionis-like objects) that also have 
a triangular H-band continuum. These objects are flagged in Table 1. As an input to this 
routine, we used the best-fit spectral type as well as the latest spectral type within the 
uncertainties. Of the 27 objects in our sample with a well matched late spectral type, 15 
(55.5% ± 14.3%) have a triangular H-band continuum. Of the 21 objects in our sample of 
the M class spectral type, 10 (47.6% ± 15.1%) have a triangular H-band continuum. These 
results show that the triangular H-band is common among Class I YSOs with late type 
spectra, yet only about half of them show this feature. Among all Class I YSOs, only about 
14% of them have a triangular H-band. Although our sample of objects are certainly all 
very young (indeed much younger than the brown dwarfs discussed above), they are likely 
to be much more massive than brown dwarfs. The objects in our sample may have evolved 
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enough (towards higher gravities) so that only about half of the late type spectra in the 
sample show this triangular H-band structure. Also, high extinction and high veiling make 
it more difficult to identify which targets have a triangular shaped H-band continuum, so 
this result should be interpreted ClS db lower limit. 



Greene & Ladal (119961 ) used the EW of Na I (2.206 fim and 2.208 /mi) plus Ca I 



(2.263 /jm and 2.266 /jm) versus the EW of CO as a gravity sensitive indicator. On this 
plot, most YSOs fall closer to the line formed by luminosity class V than luminosity class III 
stars, with FU Orionis-like stars below the line made by luminosity class III stars. Figures 4 
and 5 show the locations of our Class I protostars on this plot. The EWs for the luminosity 
class III (dwarf) and V (giant) stars were measured fro m stars in the S p eX sp ectral library 
in the same way as our YSOs. Similar to the result in iGreene &: Ladal (119961 ). most of our 
Class I YSOs also lie closer to the line formed by luminosity class V than luminosity class 
III stars (region D), with increasing veiling pushing the location of stars on this plot towards 
(0,0). Class I YSOs tend to have slightly lower gravity than dwarf stars, but not as low as 
giants. YSOs with triangular H-band continua are shown in red in Figures 4 and 5, and tend 
to be found to the upper right of the remaining YSOs. 



7. Conclusions 

We have presented the results of a spectroscopic sur vey of 110 Class I youn g stellar 



objects covering A =0.80 /im to 2.43 selected from the IConnelley et al.l (j2008al ) imaging 
survey. This atlas demonstrates the great diversity of spectroscopic properties among embed- 
ded protostars. We estimated the spectral types, extinction, veiling, and veiling temperatures 
of these protostars as the data allowed. Extinction was estimated using broadband colors 
from 2MASS, emission line ratios, and modeling of the spectrum. The median estimated 
extinction is A„ = 14.7 with a median uncertainty of 3.7 magnitudes based on continuum 
fitting, with the estimated extinctions ranging from A„ =0.0 to A„ =57.8. Scattered light 
can lead to the extinction values being underestimated. Veiling ranges from being negligible 
to completely veiling the photospheric emission. When we were able to constrain the veiling, 
the median veiling was ik = 0.48 with a median uncertainties of +1.32 and —0.24. Spectral 
types were estimated by matching photospheric line equivalent widths with values measured 
from library spectra of dwarf stars, and we present the distribution of derived spectral types. 
The median best-fit spectral type is K5 with a median uncertainty of 4 sub-classes. The best 
fit spectral types in this sample range AO to M4. Our main conclusions are: 

1) ~ 92% of embedded protostars show some sign of active accretion: Br 7 emission 
and/or high veiling. 
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2) Having divided our sample between targets with low veiling and high veiling, we 
found that targets with high veiling are over represented among targets with strong reflection 
nebulosity whereas targets with low veiling are over represented among targets with a weak 
reflection nebula or no reflection nebula. 

3) There are no embedded early type stars in the Taurus-Auriga and Perseus star forming 
regions and a relative excess of embedded early type stars in the Orion and Ophiuchus 
regions. This is consistent with what has been previously found among T Tauri stars. 

4) Observations showed that this sample includes ten FU Orionis-like objects, including 
two new embedded objects with FU Orionis-like near-IR spectra. All three FU Orionis-like 
stars with data near 1 pm show deep He I absorption. Although FU Orionis-like stars 
have high accretion luminosities, they are never observed to clearly show Br 7 emission. We 
speculate that the relationship between Br 7 line flux and mass accretion collapses at very 
high mass accretion rates. 

5) The presence of the [Fe II] and H 2 emission lines are highly correlated, suggesting 
that H 2 is also shock excited. 

6) CO was observed in emission in 27 sources. CO emission was usually only seen when 
the veiling was high, and such targets always show Br 7 emission as well. Although all FU 
Orionis-like stars showed strong CO absorption, several other stars showed CO absorption 
in excess of what is expected based on the photospheric spectral type of the star. 

7) Both components of 10 binaries were observed, with a median projected separations 
ranging from 45 AU to 4500 AU. The veiling of the components of binaries are correlated, 
with both components tending to have either low or high veiling. This is also consistent with 
what has been identified previously among T Tauri stars, but this is now observed at higher 
accretion rates and younger ages. Also, in several cases the spectra of the components of a 
binary are so similar as to be difficult to visually tell apart. 

8) Roughly half of the M-type stars in this sample have a H-band continuum with a 
triangular profile, a signature of low gravity. 
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A. Class I YSO Spectral Atlas 

Figure 13 presents the spectrum of each of the targets that we observed from 0.8 /im to 
2.45 /mi, although many objects are too faint and/or red to have any detectable flux short- 
ward of H-band. The locations of significant absorption and emission features are shown 
across the top of each page. The vertical axis is normalized by the maximum flux value in 
the spectrum of each target. Of the 110 targets that we observed, we detected the spectrum 
of all of them at K-band. We were able to detect the continuum in H-band for 107 (97%) 
targets and in J-band for 67 (61%) targets. 



B. Comments on Individual Objects 



IRAS 0322 0+3035 (N) L1448 IRS 1, RNO 13. This is a close (1"37) binary with IRAS 
03220+3035 (S) rtConnellev et al.l l2008ah . This component shows strong [Fe II] (1.644 /mi) 
emission which the spectral image shows to be spatially extended. Among other species, we 
observed emission from He I (1.083 /im, 2.058 /mi), O I (1.129 /im), [Fe II] (1.257 /mi), Fe 
II (1.688 /mi, 1.742 /im), and Na I (2.206 /im, 2.209 /im) 

IRAS 03220+3035 (S) The binary companion to IRAS 03220+3035 (N), separated 
by 1'.'37 (Connelley et al. 2007). Although this object is likely as young as any other Class 
I YSO (being the companion of a Class I YSO and also having low gravity), the spectrum 
otherwise remarkably shows few signs of youth. The spectrum is consistent with there being 
no veiling and the only emission line (Br 7) is very weak. Were this object not embedded 
and the companion of a mass accreting protostar, its spectrum could have lead to this object 
being classified T Tauri star. 



IRAS 03260+3111 (E) SVS 3. This is a close (0'.'55) binary flConnellev et al.ll2008ah . 
Light from both s tars c ontributed to this spectrum. We estimate a spectral type of F4-G0, 
Doppmann et al.l (120051 ) estimate a T e jj=5600 K, approximately G4. 
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IRAS 03260+3 111 (W) This is source 205 Per in Evans et all (120091 ) and Per-emb50 



Enoch et al.l (120091 ). which they classify as a Class I YSO. 
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IRAS 03301+3111 The best-fit spectral type is M2, but the spectrum h as a stronger 



break at 1.32 fim than a K7 star, suggesting l ow gravity. This is source 289 Per in lEvans et al. 



(120091 ) and Per-emb64 in Enoch et al.l (120091 ). which they classify as a Class I YSO 



IRAS 03445 +3242 HH 366 V LA 1 and Ba r nard 5 IRS. This is source 505 Per in 
Evans et al.l (120091 ) and Per-emb53 in Enoch et al.l (120091 ) , which they classify as a Class I 
YSO. 

IRAS 03507+3801 Source of HH 462. The CO band heads run together, and are 
not distinctly separate band heads, similar to IRAS 04073+3800. This object shows He I in 
absorption. Pa /3 is strongly seen in absorption, but Br 7 is not. 

IRAS 0401 6+2610 L1489 IRS and the source of HH 360, and a well known bipolar 
reflection nebula ( Padgett et al. 1999h. The H 2 emission is spatially extended and has been 
mapped with H 2 imaging (ILucas et al.ll2000l ). 



IRAS 04073+3800 PP 13S, source of HH 463. Thi s object has an un u sual r eflection 
nebula with a ring-like morphology (jConnelley et al.ll2007l ). ISandell fc Aspinl (119981 ) identify 
this as an FU Orionis-like star. The near-IR spectrum is dominated by a featureless contin- 
uum and deep CO absorption. However, the CO band heads for this target are unusual in 
that they do not show the usual saw-tooth pattern, but rather appear blended together and 
the band heads are not distinct. 

IRAS 04108+2803 (W) The spectrum has a stronger break at 1.32 /im than an M2 
star and a triangular H-band, suggesting low gravity. 

IRAS 04113+2758 (S) This YSO's spectrum shows a stron ger break at 1.32 11m 
than a M2 star as well as a triangular H-band, suggesting low gravity. IWhite fc Hillenbrand 
(120041 ) estimate a spectral type of M3.5+1 for this object. 



IRAS 04181+2655 Designated as IRAS 04181+2655 B in lBeckl (120071 ). who estimated 
the spectral type to be K7+1, r fc =0.0, and A„=40. Our results for the spectral type of M1^J{, 
r fc = 0.24g;j$, and A„ = 51.4^ are consistent. 



IRAS 04181+2655 (S) Designated as IRAS 04181+2655 A in iBeckl (120070 . who 
estimated the spectral type to be M3, rfc=1.3, and A t) =17.8. We derived a higher A„, a lower 
ik, and estimate that the spectral type can be from GO to M3. 



IRAS 04189+2650 (E) FS Tau A. This is a very close binary star flChen et al.lll990h . 
The H2 emission is also found to be spatially extended. 

IRAS 04239+2436 HH 300 VLA 1. The spectrum shows high veiling, with a wide 
variety of species in emission (e.g. H, H2, [Fe II], Na, Mg, CO) 



-33 - 



IRAS 04240+2559 DG Tau. Strong atomic H and Ca II emission lines. The best 
fit spectral type from photospheric lines is G3, but we also observed water absorption and 
weak CO absorption. An M2 spectrum provide s a better continuum f it in th e H-band region. 
We estimate a spectra l type of G1-K2, whereas IWhite &: Hillenbrand! (120041 ) estimate K1-K5 
and iDoppmann et al.l ( 120051 ) estimate a T e //=4000 K, approximately K7. 



IRAS 04248+2612 HH 31 IRS 2. The spectrum has stronger break at 1.32 //m than 
an M3 star and also has a triangular H-band, s u ggest ing: low gravity. We estimate a spectral 
type of M0-M4, whereas White k Hillenbrand! (J200J) estimate M5.5+1. 



IRAS 04292+2422 (E) V806 Tau Haro 6-13. We estima t e a sp ectral type of G7 (with 
acceptable fits from GO to M3), whereas White fc Hillenbrand! (12004! ) estimate MO+0.5. 



IRAS 04292+2422 (W) V806 Tau, Haro 6-13. Spectral type fitting finds a best 
match of G7 (from GO to M3), however the CO absorption is in excess of what is expected 
from a K5 dwarf star. 

IRAS 04295+2251 L1536 IRS. Has a double p eaked spectral type goo dness-of-fit 
curve: either K0-K7 or M1-M4, both need some veiling. IDoppmann et al.l (12005! ) estimate a 
T e j/=3400 K, approximately M3. Values in Table 2 are based on fits to an M2 star, which 
had the best RMS continuum fit to the data. iBecki (120071 ) estimated the spectral type to be 
K7, rfc=0.9, and A„=28. We derived similar values for r^ and A„. 

IRAS 04315+3617 This spectrum shows strong He I absorption. The calcium triplet 
near 0.85 iim is also strongly seen. 

IRAS 04530+5126 V347 Aur, RNO 33. The spectrum has a stronger break at 1.32 /xm 
than an M2 star and it also has a triangular H-band, suggesting low gravity. We also observed 
very strong blue shifted He I absorption as well as emission near rest velocity. 

IRAS 04591-0856 Gy 2-13 1. Spectral type fitting finds a best match of K7 (from K4 
to M3), however the CO absorption is in excess of what is expected from an M3 dwarf star. 
The H-band continuum is better fit with an M2 spectrum than a K7 spectrum. The spectrum 
of this object has prominent H2 emission lines, as well as broad (FWHM=640 kms -1 ) and 
red-shifted (315 kms -1 ) Br 7 emission. 

IRAS 05327—0457 (W) The best-fit spectral class is A5, and the spectral classes of 
acceptable fits range from A0 -F7. Despite this, CO bands are weakly seen in absorption. 
Adaptive optics observations (jConnelley et al.ll2008al ) have shown that this is a close (0'.'14) 
binary star, raising the possibility that the water and CO absorption features could be from 
the companion having a late-type spectrum. 
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IRAS 05404—0948 The best-fit spectral class is FO, and the spectral classes of accept- 
able fits range from A5-F3. D espite this, water and C O bands are weakly seen in absorption. 
Adaptive optics observations (jConnelley et al.ll2008al ) have shown that this is a close (0'.'16) 
binary star, raising the possibility that the water and CO absorption features could be from 
the companion having a late-type spectrum. 

IRAS 05450+0019 A new FU Orionis-lik e star, and the illuminating source of a large 
bipolar reflection nebula ( jConnelley et al.l 120071 ). 



IRAS 05513—1024 V1818 Ori. This object has an early type stellar spectrum where 
veiling dominates the spectrum at H and K-band, with CO in emission. This object also has 
strong He I absorption near rest velocity. Best fit veiling from photospheric line matching is 
r k = 8.8. 

IRAS 06297+1021 (W) This object has a spectrum that is unique in this sample. The 
spectrum of this object has several features in common with spectra from other FU Ori and 
FU Ori-like objects in this sample (e.g. high veiling, strong water absorption bands, strong 
blue shifted He I absorption). However, the spectrum also shows several strong emission 
lines, including CO and many metals (e.g. Mg & Na) that is uncharacteristic of other FU 
Ori and FU Ori-like objects in this sample. Due to these amonalies, we do not consider this 
to be a candidate FU Ori-like object. 

IRAS 06393+0913 Our photospheric line matching code found a best-fit spectral 
type of G7, with acceptable fits from GO to MO. As with other FU Orionis-like stars, the 
spectrum shows no emission lines and deep CO absorption. However, the spectrum lacks the 
deep water vapor absorption bands commonly seen in the spectra of FU Ori-like stars, and 
thus we do not classify it as such. 

IRAS 16235-2416 p Oph SI, Elias 2-25. Spectral type fitting finds a best match of 
FO (from A3 to F4), however the CO absorption is in excess of what is expected from an F4 
dwarf star. 



IRAS 16240-2430 ( E) This is also kn own as EL29, GY214, 871 Oph in Evans et al. 



( 120091 ) and Oph-embl6 in Enoch et al.l (120091 ) . which they classify as a Class I YSO. 



IRAS 16240-2430 ( W) This is a l so kno wn as WL16, GY182, 857 Oph in Evans et al. 



( 120091 ) and Oph-emb21 in Enoch et al.l ( 120091 ). which they classify as a Class I YSO. This 
object is an early type star and shows very strong CO emission. 

IRAS 16288 -2450 (E) This is also known as L1 689 IRS 5, IRS67, 1003 Oph in 
Evans et al.l ( 120091 ) and Oph-emblO in Enoch et al.l (120091 ). which they classify as a Class I 
YSO. 
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IRAS 16288— 2450W p Oph S, near the position of the source Oph-emb7 in Enoch et al 



( 120091 ). Their position is at the same RA, but is 1' to the south. 



IRAS 16289-4449 HH 57 IRS and V346 Normae. This object is considered to be 
an FU Orionis type object since it was observed to erupt. In common with other FU Ori 
and FU Ori-like objects, its spectrum has strong blue shifted He I absorption with no other 
component to the line. Unlike several other FU Ori and FU Ori-like objects, the spectrum 
of this object shows strong CO emission. The Paschen f3 line has a P Cygni profile with 
weak blue shifted absorption. 

IRAS 18250—0351 NZ Serpentis. This star shows very strong atomic hydrogen emis- 
sion, including emission from the H I Pfund series. 

IRAS 18270—0153 (6) A new FU Orionis-like s tar, and a member of a close binary 



in a small group of young stars (jConnelley et al.ll2008al ). 



IRAS 18341—0113 (N) Spectral type fitting finds a best match of M2 (with acceptable 
fits from KO to M4), however the CO absorption is in excess of what is expected from an 
M4 dwarf star. The water vapor absorption bands are much stronger than expected for a 
K5 star (the best-fit type), suggesting that the true spectral type may be later. Indeed, the 
narrow photospheric lines of a M3 type spectrum was nearly as good of a fit as the lines from 
a K5 spectrum. The H-band continuum is strongly peaked, suggesting low gravity. This is 
one of the few objects where the photospheric lines from luminosity class III stars (giants) 
are a better match to the YSO's photospheric lines than photospheric lines from luminosity 
class V stars (dwarfs). 

IRAS 18341—0113 (S) The spectrum of this source has deep water vapor absorption 
bands and lacks narrow photospheric features, in common with other FU Ori and FU Ori-like 
stars. However, the characteristicly strong CO band head absorption of FU Ori objects is 
observed to be quite weak. 

IRAS 19266+0932 Parsamian 21 and source of HH 221. Our spectral type matching 
code finds a best-fit G3 spectral class, with a range of uncertainty from F5 to K2. However, 
this source also shows water and excess CO absorption. Continuum fit results in Table 2 



are based on an M4 input spectrum, which was the best continuum fit. IStaude et al.l (119921 ) 
identified this object as an FU Orionis-like source. 



IRAS 20568+5217 HH 381 IRS, a known FU Orionis-like object (IReipurth & Aspin 



19971 ). Continuum fit results in Table 2 are based on an M8 input spectrum, which was the 



best continuum fit. 



IRAS 21352+4307 In Barnard 158. The spectral image shows that the H 2 emission 
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is spatially extended to the north by ~l'/2. 

IRAS 21454+4718 V1735 Cyg. This one of the few targets with an FU Orionis-like 
spectrum, and only one of three where we observed strong He I absorption. Skinner et al. 
(2009) found that this FU Orionis type star is the source of hard X-rays from very hot 
plasma, uncharacteristic of accretion shocks. They note that cooler plasma from accretion 
shocks could be present and heavily obscured. An embedded Spitzer sourc e 24" to the WNW, 



detected as an L'=10.7 source in images taken by lConnelley et al.l (J2008aJ), was also detected 



in X-rays. Continuum fit results in Table 2 are based on an M8 input spectrum. 



IRAS 22051+5848 HH 354 IRS, a known FU Orionis-like object flReipurth fc Aspin 



19971 ). 



IRAS 22324+4024 LkHa 233. F3 is the best-fit spectral type based on matching the 
photospheric absorption lines. Veiling based on continuum fitting is r^ = 8.4lfg, whereas 
ik = 6.2I32 is the best result from photospheric line fitting. However, the H and K band 
regions of the spectrum are best-fit by a warm (920^qq K) veiling emission. The early 
photospheric spectral type combined with the relatively cool dust suggests that there is a 
significant gap between the star and the inner edge of the dust disk. However, the presence 
of Ca II, Pa j3 and Br 7 emission shows that gas is being accreted onto the star. A solution 
to this problem is to propose that this is a close binary system with an early type primary 
star and a lower mass Class I YSO or T Tauri star with a strong IR excess as a companion. 
Early type stars are commonly in such binary systems (e.g. Baines et al. 2006, Corporon 
& Lagrange 1999). Although stars with high veiling tend to have fainter absolute K-band 
magnitudes than A or F type stars, the brighter high veiling stars have similar absolute 
K-band magnitudes as the fainter A or F-type stars. As such, it is possible for the K-band 
flux of a companion with high veiling to rival the flux from the early-type primary star. 
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Fig. 1. — The spectral index distribution of this sample of Class I protostars based on 12 /iin 
to 100 /im IRAS data. The median spectral index for the spectroscopically observed sample 
(in black) is +0.5. The red dashed histogram only includes spectroscopically observed targets 
in the Taurus, Perseus, and Auriga star forming regions. The solid gray histogram only 
includes spectroscopically observed targets in the Orion and Ophiuchus star forming r egions. 
The green dashed histogram includes all objects presented in lConnelley et al.l (j2008al ). That 
sample has a median spectral index of +0.79, suggesting that the spectroscopically observed 
sample is more evolved. 
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Fig. 2. — For Figures 2 and 3, the black histogram includes all targets, the red histogram 
only includes targets in the Taurus-Auriga and Perseus star forming regions, and the the 
green histogram only includes targets in the Orion and Ophiuchus star forming regions. This 
is an approximation of the distribution of photospheric spectral types for this sample of Class 
I YSOs. These histograms are the average the normalized goodness-of-fit curves from the 
spectral type fitting routine for all of the stars for which a photospheric spectral type could 
be estimated. 
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Fig. 3. — The distribution of the best-fit spectral types. This figure also clearly shows the 
dearth of early type Class I stars in the Taurus- Auriga and Perseus regions. 
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Fig. 4.— The equivalent width (EW) of lines of Na (2.206 and 2.209 fjm) and Ca 
(2.263 /im and 2.266 /im) versus the CO band head starting at 2.294 /im. The upper 
spectral class designations are the data for dwarf (luminosity class V) stars, and the lower 
spectral class designations are for giant (luminosity class III) stars taken from the SpeX 
spectral library. The data for the YSOs are most consistent with dwarf stars (region D) with 
some additional CO absorption, possibly due to lower gravity. Veiling tends to push the EW 
values towards zero, leading to the trail of data points towards (0,0). The location of YSOs 
which have a triangular H-band continuum, an indicator sensitive to low gravity, are shown 
in red. There are also a number of YSOs with strong CO absorption, found to the right of 
the figure (region E), which are FU Orionis-like stars. 
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Fig. 5. — This shows the full range of CO values, including those objects where CO is 
in emission (region F). When CO is seen in emission, then Na is usually also in emission. 
Calcium is not seen in emission, but in these cases the veiling is too high to see it in 
absorption. 
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Fig. 6. — The equivalent width of the [Fe II] line at 1.644 /im versus the H2 S(l) v=l-0 line 
at 2.122 fim. The red error bars are for targets with high veiling and the black error bars 
are for targets with low veiling. Targets with low veiling have a lower average [Fe II] and H 2 
EWs compared to targets with high veiling. For targets with low veiling, few targets have a 
[Fe II] EW less than -6 A and the H2 EW less than -3 A. We detected [Fe II] emission from 
all targets with H2 EW less than -3 A. 
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Fig. 7. — The equivalent width of the Br 7 line versus the H2 S(l) v=l-0 line at 2.122 //m. 
The red error bars are for targets with high veiling and the black error bars are for targets 
with low veiling. For targets with low veiling, few targets have a Br 7 EW less than -5 A and 
the H2 EW less than -3 A. The limit of the Br 7 line EW versus the H2 EW appears to be 
bound by a hyperbolic relation, where the maximum Br 7 line EW is roughly 50 divided by 
the H 2 EW. This suggests that these emission lines, and thus their excitation mechanisms, 
may be inversely correlated to each other. 
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Fig. 8. — This figure shows the relationship between the equivalent widths of the Br 7 line 
and the first CO band head at 2.293 /xm. These lines are close enough that differences 
in veiling and extinction effects are ignored. The dotted lines divide the figure into three 
regions. Region A is where the CO band heads are in emission. In this region, Br 7 is always 
seen in emission (the gray shaded area), and never in absorption. Region B is where the CO 
equivalent widths are consistent with main-sequence photospheric absorption. In this region, 
Br 7 is often seen in emission, but it can be undetected or seen in absorption. Region C is 
where the CO equivalent widths are greater than is expected from a photosphere. In this 
region, Br 7 is very rarely seen in emission, and only very weakly. Targets with high veiling 
(red) are clearly segregated from targets with low veiling (black), such that only targets with 
high veiling have CO emission or strong CO absorption (i.e. FU Orionis-like objects). 



-50 - 




r_k 

Fig. 9. — An example of the output of our continuum modeling code. The upper left panel 
shows the RMS fit contours as a function of extinction (A„) and veiling temperature. The 
upper right panel shows the same data as a surface. The lower left panel shows the best 
fit's RMS value versus the amplitude of the veiling (r^). The lower right panel shows the 
observed spectrum (black) with the modeled spectrum (green) overlaid. The red curve is a 
single temperature blackbody using the best fit veiling temperature, but not scaled according 
to the best fit r fc . This is the result for IRAS 22324+4024. 
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Fig. 10. — Histogram of our extinction estimates (A„) based on our continuum modeling 
routine (solid line) and 2MASS photometry (dashed line). 
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Fig. 11. — Plotting our estimates based on continuum modeling and 2MASS colors 
against each other. The solid line has a slope of 1 and represents the ideal location of the 
data points on this plot. The extinction estimates based on 2MASS photometry tends to be 
~ 30% lower than the extinction estimate based on continuum modeling. Error bars are not 
shown for clarity. 
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Br y Equivalent Width 



Fig. 12. — The distribution of the Br 7 equivalent widths, divided between targets with high 
veiling (black solid line) and low veiling (red dashed line). The grey solid li ne shows the Br 7 
equivalent width histogram for the T Tauri stars in iMuzerolle et al.l ( 119981 ) for comparison. 
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Fig. 13. — 0.80 to 2.43 /im spectra of our sources. 
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Fig. 13.- 



0.80 to 2.43 /zm spectra of our sources. 
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Fig. 13.- 



0.80 to 2.43 /zm spectra of our sources. 
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Fig. 13. — 0.80 to 2.43 /zm spectra of our sources. 
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Fig. 13.- 



0.80 to 2.43 /zm spectra of our sources. 
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Fig. 13. — 0.80 to 2.43 /im spectra of our sources. 
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Fig. 13.- 



0.80 to 2.43 /zm spectra of our sources. 
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Fig. 13.- 



0.80 to 2.43 /zm spectra of our sources. 
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Fig. 13. — 0.80 to 2.43 /zm spectra of our sources. 
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Fig. 13. — 0.80 to 2.43 /im spectra of our sources. 
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Fig. 13. — 0.80 to 2.43 /zm spectra of our sources. 
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Fig. 13. — 0.80 to 2.43 /zm spectra of our sources. 
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Fig. 13. — 0.80 to 2.43 /im spectra of our sources. 
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Fig. 13. — 0.80 to 2.43 /im spectra of our sources. 



-68 - 




Fig. 13.- 



0.80 to 2.43 /zm spectra of our sources. 
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Fig. 13. — 0.80 to 2.43 /zm spectra of our sources. 
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Fig. 13. — 0.80 to 2.43 /im spectra of our sources. 
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Fig. 13.- 



0.80 to 2.43 /zm spectra of our sources. 
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Fig. 13.- 



0.80 to 2.43 /zm spectra of our sources. 
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Fig. 13.- 



0.8 to 2.45 fxm spectra of our sources. 
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Fig. 13.- 



0.8 to 2.45 fxm spectra of our sources. 
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Fig. 13. — 0.8 to 2.45 /iin spectra of our sources. 
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Table 1 — Continued 
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06382+1017 


n£? /i 1 no A 

Oo 41 02. o4 


i 1 n 1 pt nn 1 

+ 10 15 02.1 


12.83 


1.00 


84.4 


1 pt f~\ „i- nnnT 

15 Oct 2007 


16 


1080 


1,3 


/ \ / ' - > ( 1 - > i /id 1 o 

06393+0913 


/ 1 / • 1 , \ no 1 o 

06 42 08.13 


i nn i n on n 

+09 10 30.0 


10.45 


1.42 


28.9 


1 t* / \ j_ nnnT 

15 Oct 2007 


347 


1800 




07018—1005(1) 


07 04 13.93 


— 10 10 13.6 


10.80 


0.35 


30.3 


17 Nov 2007 


208 


1440 




07025—1204(2) 


07 04 51.62 


— 12 09 29.9 


11.60 


1.29 


49.5 


17 Nov 2007 


138 


1920 


3 


07180—2356 


nT nn no o f* 

07 20 08.36 


n a nn no n 

— 24 02 23.0 


11.86 


0.81 


186. 


1 c "\T nnnT 

16 Nov 2007 


147 


3240 


2 


16235—2416 


1 / ■ O / ■ O .1 IT 

16 26 34.17 


c\ a no no o 

—24 23 28.3 


6.26 


1.34 


159.7 


1 a T l nnnT 

14 Jul 2007 


872 


480 


6 


1 z' 1 n /i n ri/ion/ T7^ \ 

16240— 2430(b) 


i / * ot nn /i - > 

16 27 09.43 


n A O T 1 o o 

—24 37 18.8 


7.16 


0.24 


25.6 


1 T l nnnT 

13 Jul 2007 


734 


480 


5 


-\ r> f\ a r\ n/ion/TTsr\ 

16240— 2430( W) 


i c ot nn o ■! 

lb 27 02.34 


n a ot nT n 

—24 37 27.2 


7.82 






1 T l nnnT 

13 Jul 2007 


480 


480 


5 


i coco 0,1 crn/TrA 
lDZoo— Z4oU(ii/ J 


i <^ on n i nn 

10 o2 01.00 


o A /ion 
— 24 06 4z.U 


10.13 


0.70 


5.5 


i/i T.,l onnv 
14 Jul zUO ( 


196 


1920 


1,3 


16288— 2450(W1) 


i r> 01 pt n no 

16 31 52.98 


o /I cr f o A < ■ 

—24 56 24.6 


7.85 






13 Jul 2007 


405 


480 


2,4,5? 


16288-2450(W2) 


16 31 52.98 


-24 56 24.6 


8.69 






13 Jul 2007 


334 


480 


2,4,5? 


16289-4449 


16 32 32.19 


-44 55 30.7 


7.21 


-0.04 


5.9 


13 Jul 2007 


369 


480 




16316-1540 


16 34 29.29 


-15 47 01.9 


8.28 


0.84 


11.4 


14 Jul 2007 


384 


480 


2 


16442-0930 


16 46 58.27 


-09 35 19.7 


10.92 1 


0.22 


0.7 


13 Jul 2007 


156 


1920 


2 


18250-0351 


18 27 39.53 


-03 49 52.0 




0.20 


219.7 


14 Jul 2007 


730 


40 




18270-0153(6) 


18 29 38.92 


-01 51 06.3 


11.89 






13 Jul 2007 


117 


960 


1 


18270-0153(1) 


18 29 36.69 


-01 50 59.1 


10.14 


0.49 




13 Jul 2007 


164 


1920 


2,4 


18274-0212 


18 30 01.36 


-02 10 25.6 


11.66 


0.12 




13 Jul 2007 


113 


1920 


2 


18275+0040 


18 30 06.17 


+00 42 33.6 


7.72 


-0.19 


3.4 


14 Jul 2007 


666 


480 




18278-0212(N) 


18 30 27.13 


-02 10 56.4 


11.81 






14 Jul 2007 


91 


2880 




18278-0212(S) 


18 30 27.28 


-02 11 00.2 


11.24 


1.62 




14 Jul 2007 


166 


1920 




18341-0113(N) 


18 36 46.33 


-01 10 29.5 


9.80 


0.91 




14 Jul 2007 


307 


960 


2,4 


18341-0113(S) 


18 36 46.49 


-01 10 41.6 


10.86 






13 Jul 2007 


124 


1920 


1 


18383+0059(E) 


18 40 51.87 


+01 02 12.9 


9.80 


0.50 




16 Aug 2006 


305 


960 


2 


18383+0059(W) 


18 40 51.87 


+01 02 12.9 


10.90 


0.50 




16 Aug 2006 


233 


1440 


2 


18577-3701 


19 01 08.61 


-36 57 20.1 


6. II 1 


-0.18 


1.5 


14 Jul 2007 


675 


240 


6 
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Table 1 — Continued 



IRAS a(J2000) a <5(J2000) a K s b a c L bol (L Q ) UT Date S/N d Time (s)° Flag f 



18585-3701(E) 


19 


01 


58, 


78 


-36 


57 


49 


.9 


6.60 1 


0.12 


44.3 


12 


Jul 


2007 


330 


420 




18585-3701(W) 


19 


01 


53 


68 


-36 


57 


08 


.2 


2.86 1 




44.3 


14 


Jul 


2007 


638 


80 




19266+0932 


19 


29 


00. 


.86 


+09 


38 


42 


.9 


9.68 


0.37 


3.4 


14 


Jul 


2007 


242 


960 


1 


20355+6343 


20 


36 


22. 


86 


+63 


53 


40 


.4 


10.40 


0.59 


2.5 


13 


Jul 


2007 


275 


960 




20377+5658 


20 


38 


57. 


.48 


+57 


09 


37 


.6 


9.48 


0.32 


4.8 


13 


Jul 


2007 


508 


960 


2 


20453+6746 


20 


45 


53 


.94 


+67 


57 


38 


.7 


6.47 


-0.32 


63.7 


12 


Jul 


2007 


500 


480 


2 


20568+5217 


20 


58 


21 


09 


+52 


29 


27 


.7 


7.90 


0.62 


45.6 


12 


Jul 


2007 


524 


960 


2 


21352+4307 


21 


37 


11 


39 


+43 


20 


38 


.4 


11.75 


0.17 


11.7 


14 


Jul 


2007 


294 


2400 


2 


21445+5712 


21 


46 


07. 


12 


+57 


26 


31 


.8 


10.25 


0.54 


18.5 


13 


Jul 


2007 


264 


1440 


1 


21454+4718 


21 


47 


20. 


66 


+47 


32 


03 


6 


7.03 


0.70 


106.7 


13 


Jul 


2007 


555 


240 


2 


21569+5842 


21 


58 


35 


90 


+58 


57 


22 


.8 


10.25 


0.08 


1.0 


13 


Jul 


2007 


477 


2160 


2 


22051+5848 


22 


06 


50. 


37 


+59 


02 


45 


.9 


11.23 


1.15 


73. 


13 


Jul 


2007 


88 


1920 


1 


22176+6303 


22 


19 


20 


.39 


+63 


19 


38. 


.5 


6.13 1 


0.87 


21313. 


12 


Jul 


2007 


561 


480 


1 


22266+6845 


22 


28 


02 


99 


+69 


01 


16 


.7 


10.49 


0.53 


1.8 


13 


Jul 


2007 


295 


1920 


2,3 


22272+6358(E) 


22 


28 


57. 


60 


+64 


13 


37 


.5 


8.18 


1.76 


815.5 


14 


Jul 


2007 


473 


480 


2 


F22324+4024 


22 


34 


41. 


.01 


+40 


40 


04 


.5 


9.46 


0.08 


111.6 


12 


Jul 


2007 


544 


1920 




23037+6213 


23 


05 


49. 


.76 


+62 


30 


01 


.2 


9.04 1 


1.23 


330.2 


14 


Jul 


2007 


332 


720 


2 


F23591+4748 


00 


01 


43 


25 


+48 


05 


19 


.0 


10.42 


0.60 




14 


Jul 


2007 


228 


1920 





a 2MASS coordinate for candidate YSO. 

b K-band magnitudes from Table 4 from Connelley et al. (2008), and are in the MKO photometric system unless otherwise 
noted. Magnitudes noted with 1 are from 2MASS and are in the 2MASS photometric system. 

c a is the spectral index of the source (Lada 1991), and is calculated from the IRAS fluxes of the source. 

d Peak continuum K-band signal to noise ratio. 

°Total integration time on target. 

'Flags: 1) No point source is seen at K-band. Flux likely dominated by scattered light, which affects near-IR colors; 2) 
Object associated with a reflection nebula, which may affect near-IR colors; 3) Light from both components of a known 
binary are in the spectrum; 4) Spectrum has triangular shaped H-band continuum; 5) Identified as a YSO in Table 12 of 
Evans et al. 2009; 6) Well known as a T Tauri star 
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Table 2. Spectral Type and Extinction Estimates 











Extinction (A^) 






IRAS a 


Type b 


Colors c 


Continuum 


[Fe II]° 


Paschen /3/Brackctt 7 H 2 Q(3)/S(l) 


^veiling (K) 


f 

r*= 


03220+3035(N) 




8.4 


0.6+ n ° 

w u — 0.6 


11.7 + 0.4 


18.4 + 0.1 


1220+222 

w — 140 


high 


03220+3035(S) 


M2f j 3 , 

MA 


8.4 


9.4+1° 

^ — 3.0 








o.o+2- 88 

—0.0 


F03258+3105 


M2K2 
M3 


37.2 


34 6+1'S 








12+2-2 2 

u- A — 0.0 


03260+3111(E) 


F5£ 4 

GO 


6.5 


12.6+.2'3 








0.14+2-22 

v — 0.14 


03260+3111(W) 




no O+40. 




47.82l^; 5 ^ 


680+™° 

— 180 


high 


03301+3111 


M3 


13.2 


14.7+2/ 2 , 

—3.5 


22.1 + 5.3 


15.6 ± 0.4 




0.48+°' 84 

—0.3 


03445+3242 




> 18.9 


28.8 + ?2 

^ ^ — 22 






1120+222 
u — 420 


high 


03507+3801 


G7 F5 


10.8 


J - u -°— 1.3 


11.5 + 3.9* 




1370+ 3 2° 


6+2 ,2 

u - u -0.36 


04016+2610 


5.2 


no o~l~3.2 
— 28.8 




30.711^3 


800+,, 1 " 

— 200 


high 


04073+3700 


FU Ori 


-8.4 


44.0+11 n 
— 14.0 






770+1?° 

— 370 


high* 


041 08+28031^ 


K2 K0 


20.7 


32 4+ 9 ' 6 




s.iel 3 ;^ 




12+ 1 '" 8 


04108+2803(W) 


M2&°l 

mi M4 


15.5 


17 6+I-I 








24+„ 22 


041 134-2758CN1 


M2K9 
mz M4 


14.3 


20.6±|;1 




24.7 + 0.4 


1340+222 

io1u - 640 


6+„ 2 4 
u - u — 0.0 


041 1 3+2758CS1 


M2 M0 
1V1Z M4 


10.5 


'-2.0 








oo+°- 72 

u - uu -o.oo 


04169+2702 


M2f° 


37.1 


39 2+12'?, 

— 13.0 








eo+„ 22 

0.24 


041 81 +2655("N1 


G7 G ° 


6.5 


10 n+9.5 
1J ' z -2.0 


18.0 + 4.1 


9.2 ± 0.8 


1200+ 260 
±zuu -300 


12+ 3 - 84 

u - iz -0.12 


04181+2654(M) 


Ml{S2 
ml M4 


> 50.4 


51 4+2'2 

7.5 






2000+2, 

^ vuv — 1100 


24+°-22 
w — 0.24 




C7 G0 


28.7 


43 8+ 5 ' 2 




24.69l 2 ;l 6 2 




00+ 144 

u - uu -o.oo 


04189+2650(E) 


M2K1 

mi M4 


5.4 


5 2+?'o 




5.2 ± 0.4 


1550+222 

±uau — 350 


72+Jv 22 

u * ' — 0.00 


041 89+2650fWl 




4.7 


6 0+ 6 '° 


4.6 + 1.4 


11.6 + 0.3 


1 000+ 180 
1UUU_ 220 


A Q 2 - 8 
^■ y -2.8 


04223+3700 




> 44.7 


48.0+ 4 , 
—31. 






1890ll ° 


high 


04239+2436 




20.8 


35 2+i 2 






1270+122 


hi ^h 


04240+2559 


G3Si, 


0.4 


3 l+oo 




8.5 ± 0.4 -2.26I4 J, 


1630+33° 
±uou — 500 


oo+„ 22 
u - uu — 0.00 


04248+2612 


mJ M4 


8.7 


9 2+2'2 


9.7 + 5.4 


22.0 + 1.1 




oo+2'I2 


04286+1 801 


FU Ori 


6.8 




19.3 + 0.7 






hi>h* 


04292+2422(E) 


G7 G ° 

VJI M3 


11.5 


10 9 3 T o 




14.4 + 0.5 


1610+222 

iUiU - 700 


i2+2'22 
^•- lzj — 0.00 


04292+2422(W) 


M3 


31.8 


36.8I3 3 

24 8+i 1 .; 5 








0.00+?,- 74 

— 0.00 


04295+2251 


M2f° 


22.0 






1020 + ?™ 


54 1 '2 8 ,„ 


04315+3617 




7.3 


4.8+ 2 - 2 
— 2.0 






1470+3,™ 

— 230 


3.6+2? 

— 2.1 


04325+2402 




52.2 


no q + 5-2 






2000 + ? 1cn 

wvy — 1150 


high 


04365+2535 




43.8 


30.0ll ° 






740lf 7 ° 


high 


04369+4539 




17.8 


8.0+i-g 




25.9 + 0.3 


1220l 2 «° 


high 


04381+2540 


K5 K0 


27.0 


49.21";° 






0.24l°;6« 


04530+5126 


M2 M1 
1V1Z M3 


3.4 


4-4l 2 ;? 




7.5 ± 0.5 




0-36l o ; o g 


04591-0856 




4.4 


18.512J 






1270+ 730 


n 10+0. 06 
0.18_ lg 


05155+0707 




> 57.6 


28.0+ 2 |° 






high 


05256+3049 




> 11.6 


24 0+ 18 '° 
^• u -24.0 




29.89 + 2.81 


1500l 5 7 °° 


high 


05289-0430(E) 




8.5 


4-0-4.8 


6.0 + 1.7 


16.8 + 0.7 


9201™ 


8.0l 5 3 ;° 


05302-0537 






24+|t 




33.09 + 1.90 


6301H™ 


high 


05311-0631 


M4 M2 


5.7 


42.8+^3 






18901^° 


-\ no+i.32 
± -' ; ' z -0.48 


05320-0300 


M2 K0 
mz M4 


12.4 


14.ll|1 






1650l 3 ™ 


o.i8l ;g°o 


05327-0457(W) 


AC AO 


14.6 


26.6+| : ^ 








0.48l 3 ; 3 2 « 


05357-0650 




-2.1 


6.ol : D 






1620+|° 


c: no+0.00 
o.oz_ 3 7g 


05375-0040(1) 




5.5 


6-4 +2 5 7 




5.6 ± 0.3 


1590_44Q 


00+ 1 96 
u - uu -o.oo 


05375-0040(2) 


M2^ 


7.2 


6.2t« 




5.7 + 0.4 




n 42 +0.70 
u - 4z -0.28 
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Table 2 — Continued 











T?v+i'n/>f inn / A \ 
IljXIlIlCllOIl ^.tV^ ) 








IRAS a 


Tvnp b 


VJVJ1U1 o 


(~! n ti t i ti 1 1 ii m 


We TTl e Paschen fi/Brackett -v 




T -7- d (K) 

A veiling V X V 


1 k 


05375-0040(3) 


G7 K0 


7.8 


9-OtSI 






1410lf g 


o-i4l :?2 


05378-0750(E) 




25.4 


28.0+^- 






1240+I™ 

— 4UU 


high 


05379-0758(1) 




15.2 


12.2+50 


11.4 + 3.3* 




1410111° 


5.0+2° 


05379-0758(2) 


K2 K0 


12.7 


15.6+1'| 








0.12+l ; ? 2 

— U. i-Z 


05379-0758(3) 




53.9 


+7 
32.0_22 o 










F05384-0807(l) 


M2 M4 


> 15.0 


50.4+|-| 








0.30+if 


F05384-0807(3) 


M2 M°4 


> 2.0 


36.8+I 3 4 3 








o.54i ; s 8 


F05391-0841(N) 




> 4.5 


53 2 + 9 - 






1060+2S 

— 4UU 


high 


F05391-0841(S) 




> 2.0 


4:2.0+1°- 






1120+ 8 8 80 


high 


05403-0818 




13.1 


3 2+16.0 


13.0 + 0.6 22.8 ±1.4 




1040152° 


high 


05404-0948 


FO^l 


2.5 










O.24I0I 4 


05405-0117 


M2£ 3 , 

M4 


19.3 


25.6+|? 


11.9 + 3.8 






1.56+} : ° 8 


05427-0116 


M2 M°3 


18.5 










0.36lJ°2 


05450+0019 


FU Ori 


13.5 


31 5+ 4 ' 2 








high* 


05513-1024 


f 0f° 


11.2 




11.0 + 0.4 




1420+i^ 


8.8+ai 


05555-1405(N) 




5.3 


4-Otlo 


14.0 ±0.2 




1410+Hg 


4-9+21 


05555-1405(S) 


M2 M3 


3.5 


7.8 + f 3 


-1.1 ±3.5 7.6 ±0.8 




nool 4 ™ 


o.54i ;i 4 


05598-0906(1) 




16.0 




33.1 ±0.6 




1530 cnii 
— DUU 


high 


06047-1117 




3.6 


33.6 g 






1 780+ 1030 


high 


06297+1021(E) 




15.9 


7 +i3.b 

4 


17.3 ±4.6 24.8 ±1.0 




1130+J™ 


high 


06297+1021(W) 


FU Ori 


5.0 


8.4+1* 


2.6 ±3.4* 16.4 ±1.3 




940+51° 


high* 


06382+1017 




-5.4 


5.01^6° 




24.69 ± 1.35 




high 


06393+0913 


G7 G0 

MO 


30.6 


30.2 2 f 2 








0.00+8'8o 
— u.uu 


07018-1005(1) 




11.7 




8.8 ±2.0 




1240lH° 


high 


07025-1204(2) 






2 5 . 6 n c t; 
ZO.D 






1210+32 


high 


07180-2356 




> 11.6 


4-81^8 






790li^° 


high 


16235-2416 


F0£3 


8.8 


16.2+2g 






n oc+0.06 

0.36_ 06 


16240-2430(E) 




49.3 


57.8+ 30 ' 






iiooll°° 


high 


16240-2430(W) 




29.5 




46.0 ± 0.8 






0-8181 


16288-2450(E) 




> 18.9 


30.8~^3Q"g 






840™ 


high 


16288-2450(W1) 


M2 K0 


13.5 


i7.ol 3 °6 5 


22.3 ±0.6 






r.^u_ x og 


16288-2450(W2) 


MZ M4 


13.5 


21.411-2 








n+° 84 


16289-4449 




3.9 


5.6+3H 






120011™ 


6 3+ 4 ° 


16316-1540 




2.7 


16.0tS-g 


-1.0 ±6.7* 






n oo+ 1,32 
u - uu -o.oo 


16442-0930 




8.6 


18.01^ 






12901™ 


1.411; 4 


18250-0351 




7.3 


>7 r)+3.3 
(.Z_ 7 2 


5.5 ±0.03 




1 9/tn+ 22 ° 

r^u_ 400 


Q 9 + 6.6 


18270-0153(6) 


FU Ori 




29.4l 26 5 








high* 


18270-0153(1) 




12.6 


I4.0l|;° 






" ' u -600 


i - i -0.8 


18274-0212 






14-4152:2 




23.30 ±0.93 


70O±||8° 


high 


18275+0040 




1.3 


c+3.2 
d -°-3.6 


14.9 ±0.5 




1410112° 


11. 2l 25 ' 


18278-0212(N) 






57.6l^ 7 5 6 






2000l° 500 


high 


18278-0212(S) 






40.12o. 






1500l 5 9 °° 


high 


18341-0113(N) 


1x12^° 
MZ M4 


11.6 


n 4+ 1 - 6 








Ui±z -0.00 


18341-0113(S) 


FU Ori 


12.2 


19.2l|° 








high* 


18383+0059(E) 




23.0 


90 q+4.2 
""'-11.8 






1960l 4 7 ° 


high 



Table 2 — Continued 







Extinction (A„) 




IRAS a 


Type b Colors Continuum d 


[Fe II] C Paschcn /3/Brackctt 7 H 2 Q(3)/S(l) T veiling d (K) 





18383+0059(W) 

18577-3701 

18585-3701(E) 

18585-3701(W) 

19266+0932 

20355+6343 

20377+5658 

20453+6746 

20568+5217 

21352+4307 

21445+5712 

21454+4718 

21569+5842 

22051+5848 

22176+6303 

22266+6845 

22272+6358(E) 

22324+4024 

23037+6213 

23591+4748 



G3 K2 



FU Ori 



FU Ori 



FU Ori 



K2<- 



->G0 

F3£2 
F3^> 



23.0 
1.3 

I. 06 
4.2 

-3.7 
12.1 
19.7 
19.3 
5.6 

6.8 

II. 3 
11.7 
4.7 
10.9 
21.8 
22.6 
-5.6 
13.2 
8.7 



41.0 
0.00 
9.6 



+3.5 
-32.3 
+2.1 
0.0 
+4.5 
3.9 
+2.1 



8-415.9 
-,+4.8 
-0.0 
-,+5.0 

-10.3 
-,+4.6 
J -ll. 
-1+17.0 
-4.0 

14.8±^ 

,+ 15. 

-4.2 
-1+6.0 
-7.0 
-1.5 
-1.6 
24. +l0. 



o.o: 

13.0j 

15.0^ 

4.0j 

14.8 

4.2^ 

10.0 n 

11.1 H 



24.0 

16.8; 

23.4 
38.0 

4.1 
23.8 

8.2 



18. 
+3.7 
-7. 
13. 
12.8 
■1.2 
-2.2 
+8.0 
-1.5 
+1.1 
2.1 
+7.7 
1.6 
+5.8 
2.7 



I. 2 + 4.1* 

II. + 1.6 

8.2 + 3.9 

9.3 + 2.0* 



5.5 + 0.1 





17 


1 ± 0.5 




16 


8 + 0.2 


7.0 + 3.8* 






9.6 + 4.5 


21 


5 ± 1.3 




25 


8 + 0.5 


-1.1 ± 1.9* 


17 


2 ± 0.1 



23.7 + 0.6 
32.2 ± 0.6 
8.5 + 0.5 

11.6 + 0.6 



7.27 + 0.95 



2000 
1490 
1280 
1250 



+o 

-1300 
80 
140 
170 
230 
+90 
180 



141 0+ 590 

141U_ 260 

1940l 6 8 » 
1040if 4 ° 



77n +800 
' ' u -100 
-1+200 
-170 



1050Z 



1290 



690; 



920; 



+780 
460 



-160 
-150 



-100 
-100 



1500 



+500 
750 



high 

4.8i°;| 

high 

high 
0.00 
3.9 

high 

high 

high* 

high 
5.1 

high* 

high 

high* 

high 

u - iz -0.12 

8-4i?: 2 3 

0.60+tf 
0.84+L 32 



+ 2.28 
0.00 

+30 
3.9 



T+2.6 
-3.3 



b A number in parenthesis after the IRAS designation specifies the source on the findercharts in Connelley et al. 2008a 

b Bascd on matching photospheric absorption lines. Given as best fit spectral class, then the range of well fit spectral classes as the 
uncertainty range in the superscript and subscript. "FU Ori" designates objects with a spectrum characteristic of FU Orionis-likc objects 
and/or is a known FU Orionis-like object. 

c Bascd on dereddening the 2MASS colors to the T Tauri locus as defined by Meyer et al. 1997 AJ 114 1 

d Bascd on continuum fitting with the best fit spectral type as the input (when available, otherwise M2V was used). Maximum allowed 
veiling temperature is 2000 K, which often limits the upper error bar. A veiling temperature is presented only if spectral fitting reasonably 
constrains the veiling temperature. 

e Based on the ratio of the [Fe II] 1.257 /im and 1.644 (im lines. * designates targets where the 1.644 (im line is on the edge of a comparatively 
strong absorption feature or is partially blended with the Brackett 12 emission line. 

'Based on either photospheric absorption lines or continuum fitting. "High" designated objects where we could not establish an upper limit 
on rfc because the veiling is too high. * designates FU Orionis-likc spectra where the veiling is much higher than the result of the continuum 
fitting. The dearth of photospheric lines shows that the actual extinction is likley to be much higher. 
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Table 3. Common Circumstellar Lines 



ID 


A (jim) 


Emission 


Absorption 


Comment a 


Ca II 


0.850, 0.854, 0.866 


20/33 


0/33 




He I 


1.083 


27/67 


24/67 




Fe II 


1.644 


44/107 


0/67 




H 2 


2.122 


47/110 


0/110 


1 


Br 7 


2.166 


83/110 


8/110 


2 


CO 


2.293, 2.322, 2.352, 2.383, 2.414 


25/110 


63/110 


2 



a l) The 2.122/j.m H2 line is the strongest of several H2 lines observed in the K-band; 
2) Emission is circumstellar, absorption often photospheric but can also be circumstellar 



Table 4. Emission Line Equivalent Widths 



IRAS 


Call (0.866 /tin) 


He I (1.083 fim) a 


Pa (3 (1.282 fjja) 


[Fc II] (1.644 /urn) 


H 2 (2.122 /an) 


Br 7 (2.166 ^m) 


CO (2.293 fjaaa) 


o 3 9 9 o 4- 3 n 3 ^ f tsj ^ 


Q7 ^7-1-1 n ^7 


or: /i oil oc 


Q7 QQ4-1 


i o 294-0 41 




i 7 704-0 1 Q 


Qn 01 4-1 1 ^ 


039 90-1- 30 3 ^ ( ~i 












1 71 _|_n 90 

- 1 . / izcu.zy 


074-0 fil 

y .y / ztu.oi 


P039^R4-3i o^ 












9 9c;4-0 40 


7 11 4-0 AR 


UozDU-rolll^ii/J 




n O/14-O oc 

u.y4-±u.uo 








4.U4r±U.oo 


O 7/1 _l_ O c:/1 
U. f 4±U.04 


039fi04-3i 1 1 f\AA 










-1 .OUltU.UO 


9 044-0 1 1 


1 ^94-0 44 


03301 -U31 1 1 






c: 9^4-0 90 




9 Q94-0 OS. 


9 qi _|_n OS 


^ S04-0 43 


0344^4-3949 








n 7O4-0 3fi 


9 994-0 07 


ft Q44-0 1 ^ 


i p. 9^4-0 73 


UoDU ( -f-ooUl 




1 o/14-o Qfi 
l.z4±U.oU 






O O/l 4-0 07 
-U.Z4±U.U / 


n fi/i_i_n c;7 
-U.D4±U.O ( 


1 074-0 A K 
l.U ( ±U.40 


n/im f;_i_of;i o 

U4U10-|-Z01U 










O 774-0 1 fi 
-Z.il ±\J. ID 


O Qfl-l-fl OS 




O/i n7Q_L q 2 nn 
U4U / o-|-ooUU 










n K/i 4-0 1 /i 
-U.04±U. 14 


O A 04-0 1 2 
U.4zztU.lo 


OA 7Q4-1 /IO 
ZD. 1 0±1.4U 


U41Uo-f-zoUo\li/ J 










a K7_i_n oa 
-O.O ( ±U.zD 


A n74_n oo 
-4.U f ±u.zy 


4.Do±U.DU 


U41Uo-hZoUol i VV J 












1 704-0 /I 1 
-1.1 ZitU.41 


1 1 *7A 4-0 22 
11. ( 4±U.oo 


O/l 1 T Q_l_07c;2f'M'\ 
U411o-|-Z /Do(lN ) 




i. iy±u. oo, -o.ooitu.oy 


7 Qc;4_n Qn 
- / .ooitU.oU 






Q 004-0 1 P. 
-O.ZZltU.lD 


1 2fl4-0 07 
l.oD±U.Z / 


C\A 1 1 Q_l_07c;2/'C! S \ 
U41104-Z /Oo^D j 












1 /I c;4-0 OO 
-1.40±U.zU 


1 n /in_i_o 07 
lU.4U±U.y / 


041 £0-1-9709 








o 2Q4-0 R3 
-z.oyitu.oo 


1 ^74-0 0^ 
-1.0/ ICU.UO 


1 ^04-0 00 

- 1 .ouztu.uy 


1 1 Q-l-0 9^ 

i . lyztu.zo 


U41ol-|-ZDDD(lN j 




A OQJ-O 72 
-4.UoztU. / o 






O 704-0 1 7 
-U. 1 Z±\J. 1 / 


n 204-0 of; 
-U.oyztU.UD 


O 4-0 KQ 
Z.0l±U.OO 


o/i 1 21 i ofif; t^A/H 
U41ol-f-Z0D0(lvl ) 










n /iK-i-o i o 

U.4D±U.lZ 




1 n /10-UO R7 
1U.4Z±U.D / 


O/l 121 J_Ofic;f;/'CI'l 
U41ol-rZDDD( i D ) 










1 Q7_|_n no. 

-i.o / ztu.uy 


2/11 4-0 1 A 
-0.41ltU. 14 


A 1 24-0 
4. lo±U.O0 


O/l 1 20_l_Oflc:Of'T7^ 
U41oy-hZ0DU( v ii/ ) 


-A. loztU. oo 


1 t^Q-Ln io q /iQ_i_n oi 
l.DoztU.lZ, -o.4oztU.zl 






1 Ql _|_n 1 c 
-l.olltU. ID 


1 7/14_n OA 
-1. / 4ltU.UD 


/I i Q-Uo 
4. lo±U.0U 


041 204-9(^0/' AAA 

u^ioy-t-zoou^ vv j 




nn 974-1 1 

-yu.z f in id. id 




7 ^74-0 70 


c: 1 c;-|_0 00 

-o. lomu.uy 


fl 424-0 91 
-O.^OItU.Zl 


fi34-0 fi9 
U.OOltU.OZ 


O/l O O Q _l_ Q 7fin 
U4ZZo-ro /UU 












O Qfi-LH OO 

-z.ooitu.uy 


1 1 no4-i 1 f; 
-ii.uy±i.iD 


O/l OQO_LO/1 Qfi 
U4Zoy-|-Z4oO 








-y.D4itU.D0 


1 RQ-I_n o^ 
-l.Do±U.UO 


1 Q flR4_n O/l 
-10.DDltU.Z4 


1 1 f;n4-i i q 
-ll.DUltl.lo 


O/lO/lO-LOf^^O 

U4Z4U-|-ZDoy 


-oz.ozztu.oy 


i q /i 74_n oi 
-lo.4/ ±U.ol 


1 A A Q4_n QQ 
-14.4oitU.OO 




O 7Q-UO OQ 

-U. ( oztU.Uo 


7 72-UO 07 
- 1 . ( OltU.U ( 


1 OO-UO OO 

i .zunzU.zz 


04248+2612 




-2.99±0.75 


-l.o2±0.21 




-4.89±0.17 


-2.75±0.25 


8.59±0.30 


04287+1801 








-13.77+0.61 


-1.60±0.24 


-3.32±2.04 


19.43+2.08 


04292+2422(E) 


-7.62±0.48 


0.89±0.15, -1.04±0.1 


-6.04±0.54 






-2.62±0.11 


2.31+0.31 


04292+2422(W) 












0.84±0.16 


9.02+1.18 


04295+2251 








-4.17+0.70 


-6.85±0.17 


-1.10±0.08 


5.02+1.45 


04315+3617 


-58.21±2.61 








-0.14±0.08 


-10.22±0.43 




04325+2402 








-4.55+0.69 


-1.54±0.10 


-3.01±0.12 




04365+2535 








-9.96+0.71 


-1.65±0.07 


-4.49±0.16 


3.24+1.14 


04369+2539 




2.24±0.28, -3.74±0.25 


-12.26±0.60 






-5.59±0.10 


-13.80+0.44 


04381+2540 








-13.81+3.33 


-6.65±0.49 


-2.30±0.18 


6.02+0.61 


04530+5126 


-12.42±0.18 


3.17±0.07, -0.93±0.05 


-3.89±0.14 






-2.47±0.10 


7.97+0.52 


04591-0856 










-2.12±0.09 


-1.41±0.21 


10.88+0.91 



Table 4 — Continued 







-tie i ^i.Uoo [mi) 




[rc 11J ^ 1.044 ^/m ) 


112 yz.LZZ /^mj 


p T /o ice ,,™\ 

or y ^z.ioo ^/m ) 




05155+0707 








-7.35+1.48 


-7.63+0.23 


-3.62+0.16 




05256+3049 








-8.03+0.60 


-8.66+0.23 


-6.87+0.34 


-38.83+1.51 


05289-0430(E) 


-16.58+2.87 


-6.15+1.17 


-10.52+0.94 


-3.99+0.36 


-0.98+0.14 


-2.85+0.08 




05302-0537 










-17.63+0.19 


-2.90+0.16 


-2.39+0.49 


05311-0631 












-2.79+0.10 


1.17+0.30 


05320-0300 












0.89+0.28 


7.78+0.99 


05327-0457(W) 












6.32+0.49 


1.35+0.46 


05357-0650 




0.68+0.07, -0.42+0.09 








-0.72+0.04 




05375-0040(1) 


-3.09+0.30 


3.74+0.30, -4.11+0.58 


-8.86+0.30 






-2.59+0.07 


2.43+0.17 


05375-0040(2) 


-6.64+1.16 


0.87+0.21, -2.56+0.20 


-4.91+0.27 






-1.61+0.08 


5.58+0.36 


05375-0040(3) 


-13.27+1.30 


2.38+0.35, -1.73+0.27 








-5.77+0.07 


2.09+0.42 


05378-0750(E) 








-5.80+0.92 




-3.39+0.13 


-7.51+1.04 


05379-0758(1) 










-0.47+0.12 


-6.10+0.10 


-6.67+0.83 


05379-0758(2) 












-2.35+0.19 


7.76+0.69 


05379-0758(3) 












-0.87+0.23 




05384-0808(1) 












-1.80+0.14 


5.49+0.59 


05384-0808(3) 










-0.93+0.20 


-2.73+0.15 


5.27+0.58 


05391-0841(N) 












-4.73+0.15 


3.25+0.79 


05391-0841(S) 












-3.23+0.09 




05403-0818 






-11.21+2.09 


-26.46+1.43 


-3.09+0.09 


-5.28+0.13 


-18.67+1.54 


05404-0948 












6.99+0.31 


2.28+0.71 


05405-0117 










-1.11+0.12 


-3.06+0.20 


0.94+0.91 


05427-0116 












-0.26+0.28 


8.28+1.19 


05450+0019 














29.21+1.22 


05513-1024 


-1.65+0.12 


3.62+0.10 


-4.92+0.08 






-2.02+0.10 


-5.64+0.83 


05555-1405(N) 


-13.16+2.27 


-6.07+0.94 


-15.92+0.34 




-0.61+0.10 


-8.14+0.13 




05555-1405(S) 


-6.17+0.77 


-14.56+0.81 


-3.49+0.34 


-1.77+0.44 


-0.58+0.09 


-1.09+0.08 


3.93+0.63 


05598-0906(1) 






-20.75+3.67 




-0.27+0.07 


-10.54+0.09 


-10.44+0. 71 


06047-1117 








-4.89+0.57 


-1.56+0.44 


-10.45+0.16 


-62.73+3.48 


06297+1021(E) 






-5.97+0.71 


-1.40+0.34 


-0.22+0.08 


-2.14+0.07 


-2.17+1.13 


06297+1021(W) 


-7.65+0.54 


12.37+0.23 








-3.44+0.26 


-10.65+0.50 


06382+1017 










-21.92+1.33 




1.64+6.33 


06393+0913 












0.70+0.26 


19.87+1.52 


07018-1005(1) 










-1.55+0.39 


-3.70+0.15 


-49.27+2.62 


07025-1204(2) 














1.21+0.64 



Table 4 — Continued 



TT5 A Q 
IKAo 


Uall (U.ooo ^/m ) 


rle 1 (+Uoo fim) 


ra p yL.z&z fim) 


[re iij (1.044 ^/m ) 




t~) v . _ y /o 1 aa 

rsr *y (z.loo ^m) 


( ' ( \ i<) OOQ ..w, 

KjKj (z.zyo /jjn 


07180-2356 








-34.36+2.32 


-7.01+0.21 


-1.46+0.19 


-22.60+1.37 


16235-2416 












4.62+0.21 


2.46+0.44 


16240-2430(E) 










-0.25+0.04 


-3.85+0.19 




16240-2430(W) 






-7.59+1.77 






-8.09+0.08 


-51.70+2.06 


16288-2450(E) 










-0.66+0.17 


-1.13+0.18 


2.50+0.52 


16288-2450(W1) 






-6.29+0.57 






-2.71+0.18 


1.89+0.26 


16288-2450(W2) 












-1.57+0.45 


10.90+0.77 


16289-4449 


-10.68+1.41 


14.72+0.48 






-0.91+0.14 


-2.26+0.11 


-16.76+0.68 


16316-1540 




-4.29+0.33 






-0.84+0.08 


-0.64+0.07 


6.13+0.47 


16442-0930 








-5.92+0.71 


-1.83+0.07 


-2.99+0.34 


1.69+0.74 


18250-0351 




2.90+0.24, -3.95+0.23 


-127.74+1.09 






-23.20+0.06 




18270-0153(6) 








-11.39+1.56 




0.35+0.11 


3.47+1.08 


18270-0153(1) 










-1.08+0.46 


0.58+0.23 


32.17+2.80 


18274-0212 










-50.41+0.73 


-4.35+0.25 


-17.92+1.57 


18275+0040 


-3.22+0.62 


-6.65+0.46 


-4.18+0.15 






-1.83+0.08 




18278-0212(N) 












-1.16+0.27 




18278-0212(S) 












-2.12+0.35 


0.77+0.93 


18341-0113(N) 




2.65+0.40 








-1.66+0.55 


13.54+1.48 


18341 -0113(S) 








-3.79+1.19 


-0.72+0.28 




7.03+1.59 


18383+0059(E) 












-8.11+0.16 


1.70+0.64 


18383+0059(W) 










-1.27+0.20 


-4.83+0.20 


1.33+0.60 


18577-3701 


-31.14+0.39 


1.52+0.07, -1.45+0.09 


-18.50+0.24 






-6.50+0.07 


0.73+0.48 


18585-3701(E) 




-5.60+1.02 


-11.79+0.23 






-2.27+0.17 




18585-3701(W) 


-14.39+1.28 


-5.44+0.33 


-14.06+0.29 






-2.93+0.08 




19266+0932 




3.36+0.56 








-0.18+0.28 


14.74+1.94 


20355+6343 






-3.79+0.67 


-1.06+0.21 




-1.65+0.10 


1.57+0.89 


20377+5658 






-4.70+0.24 






-2.73+0.12 


-23.32+1.27 


20453+6746 


-38.56+4.06 




-31.98+1.05 




-0.48+0.13 


-6.40+0.13 


-14.14+0.95 


20568+5217 














23.61+1.20 


21352+4307 








-6.29+0.77 


-3.70+0.08 


-4.07+0.19 


-7.59+1.00 


21445+5712 








-5.04+0.39 


-7.99+0.07 


-0.63+0.12 


0.78+0.31 


21454+4718 




10.41+0.30 








-0.94+0.34 


29.03+1.57 


21569+5842 










-1.18+0.06 


-5.12+0.08 


-7.84+0.52 


22051+5848 








-6.47+1.21 


-1.51+0.50 




36.67+1.69 


22176+6303 








-1.20+0.30 




-2.70+0.15 


-4.36+0.77 



Table 4 — Continued 



IRAS 


Call (0.866 /urn) 


He I (1.083 ^m) a 


Pa & (1.282 pi) 


[Fe II] (1.644 ,um) 


H 2 (2.122 pi) 


Br 7 (2.166 fim) 


CO (2.293 pi) 


22266+6845 






-5.06+0.60 


-1.77+0.38 


-2.15+0.09 


-1.82+0.08 


2.72+0.85 


22272+6358(E) 






-8.76+0.77 






-2.42+0.09 


-5.04+0.68 


F22324+4024 


-6.59+. 21 


3.01+0.14, -0.92+0.14 




-1.64+0.18 


-0.42+0.08 


-2.19+0.06 




23037+6213 










-0.32+0.08 


-2.90+0.09 




F23591+4748 




-6.95+1.36 


-6.57+0.43 




-0.70+0.14 


-2.81+0.19 


4.58+0.48 



a The He I line is often seen in emission and absorption simultaneously, in which case positive and negative equivalent width values are given for the absorption 
and emission components. 
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Table 5. Nebulosity and Spectral Category 



# % a 

Group 1: Strong Nebulosity 

Low veiling 1/15 6.7%+l 3 f^ 

High veiling 8/15 53.3%+^;^° 

FU Ori 6/15 40.0%+^;^ 

Group 2: Some Nebulosity 

Low veiling 19/41 46.3%+|"^ 

High veiling 19/41 46.3%+|;^° 

FU Ori 3/41 7-3%±f$% 

Group 3: No Nebulosity 

Low veiling 29/54 53.7%+^° 

High veiling 24/54 44.4%+^° 

FU Ori 1/54 1.9%+*;^°° 



a l a uncertainties are presented, derived us- 
ing binomial statistics with the calculator at 
http: / / statpages.org/confint.html 
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Table 6. FU Ori and FU Ori-like Objects, and Objects with Excess CO Absorption 



IRAS CO (A) a L bot (L ) Notes 



04073+3800 


26 


.73 


22.6 


known FU Ori 






04108+2803W 


11 


.74 


0.7 


known YSO, no association as an FU Ori 


04286+1801W 


19 


.43 


20.2 


known FU Ori 






04292+2422W 


9 


.02 




excess CO absorption 


for 


G0-M3 star 


04591-0856 


10 


.87 


0.90 


excess CO absorption 


for 


K4-M3 type 


05327-0457W 


1. 


.19 




excess CO absorption 


for 


A0-F7 star 


05404-0948 


2 


.29 




excess CO absorption 


for 


A5-F3 star 


05450+0019 


29 


.20 


27.6 


new FU Ori 






06393+0913 


19 


.86 


8.9 


excess CO absorption 


for 


G0-M0 star 


16235-2416 


2 


.46 




excess CO absorption 


for 


A3-F4 star 


16289-4449 


-17 


.45 


5.9 


known FU Ori 






18270-0153W 


32 


.17 




new FU Ori 






18341-0113N 


13 


.54 




excess CO absorption 


for 


K0-M4 type 


18341-0113S 








new FU Ori? 






19266+0932 


14 


.74 


3.4 


known FU Ori 






20568+5217 


23 


.61 


45.6 


known FU Ori 






21454+4718 


29 


.03 


106.7 


known FU Ori 






22051+5848 


36 


.67 


73.0 


known FU Ori 







a CO equivalent width in Angstroms 



-89 - 



Table 7. Binary Properties 



IRAS 


D(pc) 


r(") 


r(AU) a 


AK 


Veiling 


03220+3035 


290 


1.37 


397 


0.26 


mixed 


F04108+2804 


140 


21.64 


3029 


0.69 


low 


F04113+2758 


140 


3.97 


555 


0.07 


low 


F04189+2650 


140 


19.75 


2765 


3.75 


mixed 


05375-0040 


470 


6.40 


3008 


1.77 


low 










0.94 


low 


F05391-0841 


480 


9.49 


4555 


0.05 


high 


16288-2450 


160 


2.98 


476 


0.84 


low 


18278-0212 


600 


4.47 


2682 


0.57 


high 


18383+0059 




2.15 




1.11 


high 



a Projected separation on the sky 



